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ABSTRACT 

A  computer-controlled  data  aquisition  and  analysis 
system  was  designed,  built,  and  tested.   All  of  the  required 
elements  of  a  complete  data  aquisition  system  were 
assembled,  including  transducers,  amplifiers,  signal 
conditioners,  display  devices,  computer  interface,  computer, - 
and  related  software.   A  digital  multiplexer  was  constructed 
to  multiplex  up  to  16,  16-bit  channels  for  input  into  the 
computer.   A  computer  program  was  written  to  bring  the 
various  parts  of  the  system  together.   An  experiment. 
Analysis  of  an  Airfoil  by  Pressure  Distribution,  was  used  as 
a  vehicle  to  test  the  system.   The  final  output  of  the  data 
aquisition  system,  including  graphical  information,  compared 
favorably  with  previous  results  from  an  older  data 
aquisition  system  in  use. 
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I.  INTRODUCTION 

A.   GENERAL  DISCUSSION 

The  object  of  any  data  aquisition  and  analysis  system  is 
to  collect  the  data,  process  them  in  a  desired  fashion,  and 
record  the  results  in  a  form  suitable  for  storage, 
presentation,  or  additional  subsequent  processing.   The 
simplest  method  of  data  aquisition  might  consist  of  a 
person,  reading  an  instrument,  and  recording  the  observation 
on  a  data  sheet.   Complex  data  aquisition  systems  could 
incorporate  the  aquisition  of  thousands  of  data  points, 
coupled  with  real-time  analysis  and  final  output  on  the 
device  of  your  choice. 

Figure  1  illustrates  the  basic  building  blocks  to  most 
data  aquisition  systems.   Data  aquisition  begins  with  a 
physical  property  to  be  measured,  such  as  mass,  velocity, 
pressure,  position,  etc.   A  transducer  is  used  to  convert 
the  measurable  quantity  into  an  electrical  signal  which  is 
used  as  a  direct  representative  of  that  physical  property. 
The  output  of  the  transducer  could  be  in  the  form  of  a 
varying  voltage,  current,  or  frequency.   From  this  point 
onward  until  output  is  desired  (since  human  beings  cannot 
"read"  electrical  signals  directly) ,  the  value  of  the 
original  physical  property  remains  in  electrical  form.   An 
input  circuit  is  used  to  ensure  transducer  compatibility  to 
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Figure  1 


the  rest  of  the  electrical  circuitry.   Signal  conditioning 
refers  to  various  techniques  used  to  reduce  the  effects  of 
noise.   Transmission  involves  transporting  the 
representative  electrical  signal  from  one  location  to 
another.   It  really  could  be  located  at  any  point  in  the 
block  diagram  since  all  signals  travel  an  arbitrary  distance 
between  processing  functions.   Processing  involves  molding 
the  raw  signal  into  a  form  that  could  be  used  to  make 
comparisons  with  other  signals  of  the  same  kind,  or  signals 
coming  from  different  sources.   Processing  also  includes 
data  manipulation  such  as  mathematical  calculations,  and 
type  conversions  such  as  analog-to-digital  or  current-to- 
voltage.   Amplification  could  be  used  in  the  input  circuit, 
signal  conditioning,  transmission,  or  processing  steps. 
Lastly,  the  worked  data  (or  if  desired,  the  raw  data)  is 
stored  or  displayed  so  that  it  could  be  interpreted 
immediately  through  paper  or  crt  display  output,  or  used  at 
a  later  time. 

B.   COMPUTER  ASSISTED  DATA  AQUISITION 

The  digital  computer  presents  some  unique  opportunities 
and  techniques  for  data  aquisition  and  analysis,  and  its  use 
incurs  many  benefits.   First,  it  has  multipurpose  and 
multitasking  capabilities.   It  can  record,  manipulate  or 
process,  display  and  store  data  in  a  very  efficient  manner. 
It  not  only  can  receive  information,  but  it  can  distribute 
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information  to  control  the  various  processes  involved  from 
simple  on-off  operations  to  complicated  robotics.   Second, 
it  can  make  decisions,  outlined  in  software,  that  allow  it 
to  be  effective  over  a  wide  range  of  input  parameters. 
Third,  it  is  fast.   Although  the  term  "fast"  is  relative, 
and  "fast"  is  getting  faster  every  day,  in  comparison  to  the 
capabilities  of  a  person  to  perform  the  same  duties,  it  has 
great  speed.   Fourth,  it  excels  at  repetitive,  boring  tasks 
that  often  compromise  the  routine  of  an  experiment,  without 
forgetting  steps.   Fifth,  it  is  relatively  inexpensive  in 
comparison  with  most  high  grade  experimentation  equipment, 
especially  considering  what  it  accomplishes.   Lastly,  it 
frees  its  human  operator  from  data  collecting  to  data 
interpretation,  which  is  the  task  people  (vs  computers) 
excell  at. 

Using  a  computer,  however,  does  not  solve  every  problem. 
First,  no  computer  ever  designed  an  experiment.   A  human 
must  perceive  the  need  for  information,  determine  what  he 
needs  to  know,  and  how  he  is  to  obtain  that  information.   A 
person  must  build  the  separate  components  to  a  data 
aquisition  system  and  ensure  their  compatibility.   A  person 
is  not  relieved  of  the  basic  calibration  and  other  data 
processing  fine  tuning  that  ensures  correct  results. 
Second,  "Garbage  in  -  Garbage  out."  is  always  true  with 
computers.   Errors  can  multiply  quickly,  or  worse,  can  go 
unnoticed  if  the  the  results  are  within  expectations. 
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Third,  since  transistor  logic  is  standardized,  the  data  must 
always  be  converted  into  a  form  that  computer  circuits  can 
recognize.   Fourth,  the  computer  cannot  interpret  data 
beyond  software  algorithms.   That  is,  it  cannot  decide  if 
the  results  are  meaningful.   In  the  final  analysis,  no 
computer  can  guarantee  meaningful  results,  it  can  only  make 
it  less  tedious  to  obtain  them. 

C.   THESIS  GOALS 

It  was  the  goal  of  this  thesis  project  to  design,  build, 
integrate  the  required  system  components,  and  operate  a 
complete  data  aquisition  system.   The  goal  also  included  the 
writing  of  a  computer  program  to  aquire,  analyze  and  display 
the  data  on-line,  and  store  the  required  information  for 
later  processing.   Additionally,  the  computer  program  must 
produce  a  hardcopy  report  to  include  tabular  data  and 
graphical  output.   A  design  requirement  was  that  the  system 
be  as  general  as  possible  so  that  it  could  be  used  for  a 
variety  of  tasks. 

As  a  test  case,  a  pre-designed  experiment:  Airfoil 
Performance  By  Pressure  Distribution  [Ref .  1] ,  was  chosen  as 
a  vehicle  to  test  the  system.   Measurement  requirements 
included  pressure,  temperature,  angle  of  incidence.   Data 
processing  requirements  included  input  signal  conditioning, 
analog-to-digital  conversion,  and  mathematical  data 
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reduction.   Final  output  requirements  included  on-line 
displaying  of  raw  data,  reduced  data  and  graphed  data,  and 
X-Y  plotting  of  calculated  results  with  paper  hardcopy. 
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II.   WIND  TUNNEL  AND  INSTRUMENTATION 

A.   WIND  TUNNEL  GENERAL  DESCRIPTION 

A  complete  description  of  the  Department  of  Aeronautics 
Low-Speed  Wind  Tunnel  Facility  is  found  in  Reference  1. 
Figure  2  shows  the  general  layout  of  the  tunnel.   The  wind 
tunnel  has  two  static  pressure  manifolds  that  measure  the 
static  pressure  in  the  settling  chamber  and  the  static 
pressure  in  the  test  section.   The  manifold  in  the  settling 
chamber  is  located  just  aft  of  the  screens  [Fig.  2].   The 
other  manifold  is  located  at  the  entrance  to  the  test 
section.    Pressure  from  these  manifolds  are  used  to  provide 
the  total  pressure  (with  appropriate  corrections)  and  the 
free  stream  static  pressure. 

The  airfoil  tested  was  a  NACA  66 (215) -216  section. 
Figures  3  and  4  show  the  general  shape  of  the  airfoil  as 
well  as  the  locations  of  the  static  ports.   The  3  6  static 
ports  are  tapped  from  leading  edge  to  trailing  edge,  on  both 
the  upper  and  lower  surfaces. 

B.   INSTRUMENTATION 
1.   Scanivalve 

One  48-port  scanivalve  was  used  to  measure  each  of 
the  36  airfoil  pressures,  the  settling  chamber  static 
pressure,  and  the  free  stream  static  pressure.   The 
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scanivalve  was  used  because  it  provides  an  economical  means 
of  sensing  the  pressure  of  a  multiple  number  of  sources.   It 
also  simplifies  calibration  procedures  in  that  only  one 
transducer  needs  to  be  calibrated. 

The  scanivalve  mechanism  puts  out  a  7-bit  binary- 
coded  decimal  (BCD)  signal  that  corresponds  to  the  port  (1 
of  48)  currently  connected  to  the  scanivalve  transducer. 
This  arrangement  allows  remote,  electronic  monitoring  of  the 
configuration  of  the  port  assembly  of  the  scanivalve.   This 
BCD  signal  is  logic-reversed.   That  is,  a  "1"  put  out  by  the 
scanivalve  would  correspond  to  a  "0"  in  TTL  logic.   For  this 
reason,  the  scanivalve  BCD  output  is  inverted,  bit  by  bit, 
before  it  is  used  for  input  into  Hewlett-Packard  digital 
indicators.   The  digital  indicators  provide  a  visual 
indication  of  the  configuration  of  the  scanivalve,  apart 
from  the  automatic  data  aquisition  process. 
2 .   Settling  Chamber  Static  Pressure 

Settling  chamber  static  pressure  was  taken  from  the 
static  pressure  manifold  located  in  the  settling  chamber. 
The  manifold  consists  of  four  static  pressure  ports  located 
on  the  floor,  ceiling,  and  side  walls  of  the  tunnel.   Figure 
5  shows  the  general  design  of  the  manifold.   All  four  static 
pressure  port  are  interconnected  to  reduce  the  effects  of 
flow  peculiarities. 
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3 .  Free  Stream  Static  Pressure 

Tunnel  free  stream  static  pressure  was  taken  from 
the  other  pressure  manifold  located  at  the  entrance  to  the 
test  section.   It  was  constructed  in  the  same  manner  as  the 
manifold  for  the  settling  chamber  static  pressure,  its 
location  being  the  only  difference. 

4 .  Tunnel  Dynamic  Pressure 

Tunnel  Q  indication  was  provided  by  a  manometer, 
reading  in  centimeters  of  water.   The  manometer  has,  for  its 
inputs,  the  pressures  from  the  two  static  pressure  manifolds 
described  earlier.   The  manometer  was  a  part  of  the 
established  wind  tunnel  facility. 

5 .  Temperature  Sensing 
Copper-constantan  thermocouple  probes  were 

manufactured  to  provide  inputs  for  the  tunnel  total 
temperature,  and  the  ambient  temperature.   The  probe 
measuring  the  total  temperature  was  mounted  in  the  tunnel 
wall  of  the  settling  chamber,  just  above  the  entrance  door. 
The  probe  measuring  the  ambient  temperature  was  mounted  on  a 
shelf  built  into  the  exterior  of  the  tunnel  wall  between  the 
settling  chamber  and  the  test  section.   Copper-constantan 
thermocouples  were  chosen  for  their  sensitivity  in  the 
anticipated  temperature  ranges  of  tunnel  operation. 

6.  Digital  Thermometer 

A  digital  thermometer  was  used  to  provide  an  output, 
in  degrees  fahrenheit,  from  the  two  copper-constantan 
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thermocouples.   This  panel-type  thermometer  was  mounted  in 
the  instrumentation  rack  and  had  facilities  for  manual 
operation  as  well  as  computer-controlled  operation.   It  was 
desirable  to  have  a  visible  means  of  temperature  indication, 
independent  from  the  automatic  data  aquisition  process,  so 
that  an  experimenter  could  check  the  system  while  the 
computer  was  involved  in  other  tasks. 

7.  Calibration  Pressure 

Calibration  pressure  was  provided  by  a  manually 
operated  u-tube  manometer.   The  manometer  had  the  capability 
of  providing  +-  60  cm  of  water  of  pressure  for  calibration. 
Cost  considerations,  lack  of  time,  and  the  ready 
availability  of  the  u-tube  manometer  prompted  the  choice  to 
use  it. 

8 .  Calibration  Pressure  Measurement 

Calibration  pressure  was  measured  by  a  differential 
pressure  transducer.   Its  high  input  was  connected  to  the 
calibration  manometer,  and  its  low  input  was  vented  to 
atmospheric  pressure.   This  transducer  was  also  a  part  of 
the  wind  tunnel  facility  and  was  available  immediately. 

9 .  Ambient  Pressure 

Ambient  pressure  was  measured  by  an  absolute 
pressure  transducer  that  was  simply  vented  to  the  ambient 
pressure  of  the  room  in  which  the  wind  tunnel  was  located. 
(Note:   It  was  the  intent  of  the  design  process  for  aquiring 
the  various  component  transducers  to  use  some  very  accurate 
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Digiquartz  pressure  transducers  coupled  to  a  pressure 
computer  manufactured  by  Paroscientif ic.   The  delay  of  the 
delivery  of  these  precise  instruments  prompted  the  use  of 
the  above  mentioned  transducers  to  measure  calibration 
pressure  and  ambient  pressure.) 

10 .  Angle  of  Incidence 

The  angle  of  incidence  was  measured  by  a  ten-turn 
potentiometer,  coupled  to  the  drive  system  of  the  wind 
tunnel  turntable.   Figure  6  shows  the  general  layout  of  the 
turntable  system.   The  wind  tunnel  turntable  consisted  of  a 
DC-motor  driving  a  large  sprocket  via  chain  drive.   A  shaft, 
concentrically  attached  to  the  large  sprocket  provided  the 
means  to  turn  the  wing  model.   A  simple  relay  system 
operated  the  DC-motor.   The  potentiometer  was  attached  to  a 
gear  connected  to  the  turntable.   Its  output  was  in 
millivolts. 

11.  Computer  System 

An  IBM  PC-AT  was  provided  to  accomplish  the 
necessary  data  aquisition  and  control  tasks.   A  device 
called  the  "Data  Aquisition  And  Control  Adapter"  (henceforth 
called  the  "adapter")  was  provided  to  interface  the  computer 
to  the  analog/binary  "outside  world".   The  adapter  consisted 
of  four  analog  input  channels,  two  analog  output  channels, 
one  16-bit  input  binary  channel,  one  16-bit  output  binary 
channel,  a  delay  device,  and  a  counter.   Two  adapters  were 
installed  inside  the  computer.   They  are  labelled  adapter 
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#0,  and  adapter  #1.   Analog  input  requirements  were  +  -  10 
volts  in  the  analog  input  channels.   Binary  input 
requirements  were  standard  TTL  logic. 
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III.   COMPUTER  -  INSTRUMENT  INTEGRATION 

A.  DIGITAL  PANEL  METERS 

The  design  goal  of  this  thesis  was  to  use  the 
Paroscientif ic  digiquartz  transducers  and  the  Paroscientif ic 
pressure  computer.   The  interface  of  this  computer  to  the 
adapter  would  have  been  binary.   For  a  temporary  measure, 
generic  digital  panel  meters  (DPM)  were  utilized  to  function 
as  substitutes.   The  Analog  Devices  DPM  model  #AD2024  was 
chosen  for  this  task.   It  had  a  parallel  BCD  interface  that 
was  similar  to  that  of  the  Paroscientif ic  unit.   These  DPM's 
had  slow  conversion  rates,  but  since  they  were  installed  to 
provide  indications  for  the  calibration  and  ambient 
pressure,  and  since  this  information  was  not  time  critical, 
the  DPM's  were  considered  appropriate  for  the  task.   A  DPM 
was  also  used  to  provide  an  indication  of  the  angle  of 
incidence. 

B.  MULTIPLEXER  DESIGN  AND  CONTROL 

BCD  encoded  information  is  transported  in  parallel 
fashion  (vs  serial  communications) .   Each  digit  requires 
four  bits.   Four  full  digits  require  16  simultaneous  binary 
signals  to  convey  the  information.   Since  the  two  data 
aquisition  and  control  adapters  contain  two  16  bit  channels 
between  them,  it  quickly  became  apparent  that  32  bits  would 
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not  be  sufficient  to  contain  all  the  BCD  channels  required 
from  the  four  DPM's,  and  the  scanivalve  port  location 
indicators.   A  digital  multiplexer  was  designed  to  multiplex 
up  to  16,  16  bit  channels  into  one  16  bit  input  into  the 
computer. 

The  Intel  8212  8-bit  input/output  port  was  chosen  as  the 
prime  component  of  the  multiplexer.   It  allows  a  single 
digital  input  to  control  the  throughput  of  data  signals  from 
its  inputs  to  its  outputs.   For  example,  if  a  "1"  is  applied 
to  its  control  input,  data  flows  from  its  eight  data  inputs 
to  the  eight  data  outputs.   Figure  7  illustrates  this.   If  a 
"0"  is  applied  to  the  control  input  of  the  8212,  then  the 
data  outputs  "float"  regardless  of  the  signals  at  the  data 
inputs. 

Two  8212  chips  complete  a  16  bit  input  channel  into  the 
multiplexer.   They  are  controlled  as  a  pair.   Overall 
control  is  maintained  as  follows:  select  one  16  bit  channel 
by  sending  a  "1"  to  the  pair  of  8212  chips  controlling  that 
channel,  and  by  sending  a  "0"  to  all  the  other  8212 's 
thereby  preventing  all  other  channels  from  interfering  with 
the  selected  channel.   Information  can  then  pass 
unambiguously  from  any  one  16  bit  channel  into  the  16  bit 
digital  input  of  one  of  the  computer's  data  aquisition 
adapters  (the  other  adapter's  16  bit  digital  input  channel 
was  not  used) . 
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Overall  control  of  the  mulitplexer  is  performed  through 
one  of  the  computer  data  aquisition  adapter  16  bit  digital 
outputs.   Thus  sending  a  "2"  will  enable  the  second  channel 
to  send  information;  sending  an  "8"  will  enable  the  fourth 
channel  to  send  information,  and  so  forth.   This  digital 
multiplexer  then  is  the  key  controlling  device  to  obtain 
information  from  the  DPM's,  and  the  scanivalve  port 
indicators. 

C.   SIGNAL  PATHS 

Figure  8  shows  the  general  layout  for  signal  paths  from 
each  transducer  to  the  computer.   There  are  four  different 
signal  paths  that  characterize  the  figure,  the  path  for 
scanivalve  output,  the  path  for  scanivalve  port  location, 
the  path  for  calibration  pressure,  ambient  pressure,  and 
angle  of  incidence,  and  the  path  for  temperature  signals. 

1.   Scanivalve  Signal  Path 

The  signal  path  for  the  scanivalve,  containing 
information  about  wing  port  pressure,  free  stream  static 
pressure,  and  settling  chamber  static  pressure,  passes 
through  the  signal  conditioner/amplifier.   The  signal  is 
conditioned  by  a  low  pass  filter  with  a  cutoff  of  10  hz . 
This  removes  most  of  the  noise,  including  power  line 
interference.   The  amplifier  amplifies  the  signal  from  the 


28 


^ 

Analog 
Input 

(D 

13 

Cl 

E 
o 

O  - 

Binary 
Input 

Scanivalve 

► 

Amplifier 

^ 

p 

0) 

X 

d 

— ► 

Scanivalve 
Port  Location 

^ 
w 

Calibration 

Pressure 

Transducer 

► 

Amplifier 

-^^ 

DPM 

Atmospheric 

Pressure 

Transducer 

^ 

Amplifier 

w 

DPM 

^ 

w 

Tunnel  Total 

Temperature 

Thermocouple 

DPM 

r 

Ambient 

Temperature 

Thermocouple 

Turntable 
Potentiometer 

► 

Amplifier 

^ 

DPM 

F 

w 

Signal  Paths 
Figure  8 

29 


millivolt  range  of  the  output  of  the  transducer  to  the  +  10 
volt  range  of  the  adapter.   The  conditioned,  amplified 
signal  is  input  into  analog  input  channel  #0  of  adapter  #0. 

2 .  Signal  Path  Through  a  PPM 

Millivolt  signals  from  the  calibration  transducer, 
ambient  pressure  transducer,  and  angle  of  incidence 
potentiometer  are  fed  into  the  signal  conditioner/amplifier. 
The  gain  of  each  is  adjusted  to  calibrate  each  corresponding 
DPM  in  the  desired  engineering  units.   For  example,  the  gain 
for  the  angle  of  incidence  signal  was  adjusted  at  the 
amplifier  so  that  the  angle  of  incidence  DPM  reads  out 
directly  in  degrees.   The  gain  of  the  calibration  signal  was 
adjusted  to  allow  its  DPM  to  read  in  centimeters  of  water. 
Lastly,  the  gain  of  the  ambient  pressure  signal  was  adjusted 
so  that  its  DPM  reads  in  lbs/in^.   All  three  signals  are 
conditioned  through  the  low  pass  filter  with  a  cutoff  of  10 
hz.   The  DPM  for  each  signal  has  a  TTL  compatible  parallel 
BCD  output.   The  computer  receives  these  BCD  signals  via  the 
multiplexer  into  the  16-bit  binary  input  of  adapter  #0. 

3 .  Temperature  Signal  Path 

The  outputs  from  the  thermocouples  were  input 
directly  into  a  special  DPM,  the  Analog  Devices  AD2036 
Digital  Thermometer.   This  DPM  contains  the  amplification, 
linearization,  cold  junction  compensation  necessary  for 
accurate  temperature  measurements  from  the  copper-constantan 
thermocouples.   The  DPM  outputs  directly  in  degrees 
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fahrenheit  and  also  has  a  TTL  compatible  parallel  BCD 
output.   The  binary  output  from  the  AD2  03  6  is  input  into  the 
16-bit  binary  input  of  adapter  #0,  via  the  multiplexer. 

4 .  Port  Location  Signal  Path 

The  computer  receives  scanivalve  port  location 
information  from  the  same  bit-inverters  that  generate  the 
correct  signal  for  the  Hewlett-Packard  port  indicators. 
This  BCD  signal  path  also  includes  the  multiplexer.   The 
signal  comes  into  the  computer  via  the  16-bit  binary  input 
port  of  adapter  #0. 

5.  Instrumentation  Rack 

Figure  9  is  a  diagram  of  the  instrumentation  rack. 
All  interconnections  are  made  from  the  rear  of  the  rack. 
Amplifier  adjustments  are  made  from  the  front  of  the 
amplifiers. 

D.   SYSTEM  CONTROL 

1.   Scanivalve  Control 

Figure  10  shows  the  general  control  paths  for  the 
scanivalve  and  turntable  position.   The  scanivalve  consists 
of  two  separate  units:  one  unit  contains  the  transducer, 
motor  drive,  and  port  assembly,  the  second  unit  houses  the 
solenoid  controller.   Control  of  the  scanivalve  is 
accomplished  through  the  solenoid  controller.   The  solenoid 


31 


MULTIPLEXER 


1     1    1     1 
1     1   1     1 
I     1   1     1 

CALIBRATION 


+  55.00 


PRESSURE 


ATMOSPHERIC 


+  14.70 


PRESSURE 


ANGLE 
OF 


PORT 


TEMPERATURE 


+002.0 


+  67.1 


o 


48 


INCIDENCE 


SCANIVALV 
CONTROLLER 


SCANT VALVE 
RELAYS 


AMPLIFIERS/ 
SIGNAL  CONDITIONERS 


Instrumentation  Rack 


Figure  9 


32 


SCAN I VALVE 


t 


SCANIVALVE  CONTROL 

STEP 


SOLENOID 
CONTROLLER 


STEP-UP 
RELAYS 


t 


HOME 


7 -BIT  PORT 


LOCATION  SIGNAL 


MULTIPLEXER 


BINARY 
OUTPUT 


BINARY 
INPUT 


COMPUTER 


ANGLE  OF  INCIDENCE 
CONTROL 


INCREASE  AOA  —/ 

^ 

BINARY 
OUTPUT 

DC  MOTOR 

^ 

DC-MOTOR 
CONTROL  RELAYS 

4                        ► 

^ 
^ 

^ 
^ 

^ 

^             ► 

▼ 

DECREASE  AOA  ^ 

TURNTABLE 
POTENTIOMETER 

COMPUTER 

^ 

r 

' 

BINARY 
INPUT 

^ 

DPM 

k. 

MULTIPLEXER 

k 

AMPLIFIER 

w 

r 

r 

System  Control 
Figure  10 


33 


controller  allows  two  commands,  STEP  and  HOME.   STEP  will 
"step"  the  scanivalve  one  port  location.   The  HOME  command 
sends  the  scanivalve  to  port  number  48. 

Computer  control  is  accomplished  through  the 
solenoid  controller.   It  will  accept  a  15  volt  pulse  to 
either  STEP  or  HOME.   A  pulse  is  sent  by  using  one  of  the 
16-bit  digital  outputs.   Bit  position  one  (or  the  least 
significant  bit)  of  the  16-bit  binary  output  port  of  adapter 
#0  is  reserved  for  the  STEP  command.   STEP  is  accomplished 
by  outputting  a  binary  "1",  followed  by  a  binary  "0". 
Similarly,  bit  position  two  is  reserved  for  the  HOME 
command.   HOME  is  effected  by  outputting  a  binary  "2", 
followed  by  a  binary  "0".   A  binary  "2"  will  make  the  second 
bit  position  go  high,  while  leaving  all  others  low.   Since, 
for  the  computer,  a  binary  "1"  or  "high"  is  about  five 
volts,  an  intermediate  relay  device  was  built  to  step  up  the 
voltage  (to  15  volts)  into  the  solenoid  controller. 
2 .   Angle  of  Incidence  Control 

The  angle  of  incidence  is  controlled  through  the 
turntable  assembly  located  beneath  the  test  section.   The  DC 
motor  is  controlled  through  a  relay  system  that  allows 
manual  control  as  well  as  computer  control.   Figure  10  shows 
the  general  arrangement. 

Bit  position  three  of  the  16-bit  binary  output  port 
of  adapter  #0  was  reserved  for  angle  of  incidence  - 
increase.   By  outputting  a  binary  "4",  the  relay  system  will 
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drive  the  DC  motor  in  the  positive  direction.   Similarly, 
bit  position  four  is  reserved  for  angle  of  incidence  - 
decrease.   By  outputting  a  binary  "8",  the  relay  system  will 
drive  the  DC  motor  in  the  negative  direction.   A  binary  "0" 
will  cancel  any  previous  commands  and  stop  the  turntable 
drive  motor. 

3 .   Temperature  PPM  Control  - 

The  AD2  03  6  digital  thermometer  was  different  from 
the  other  DPM's  in  that  it  allowed  six  different  inputs  as 
opposed  to  one  for  the  other  DPM's.   The  temperature  DPM 
also  has  provisions  for  manual  control  as  well  as  for 
computer  control. 

Each  temperature  channel  was  selectable  through  a 
BCD  input.   A  binary  "0"  selected  channel  0,  a  binary  "5" 
selected  channel  5,  and  so  forth.   On  the  binary  output  port 
of  adapter  #0,  bit  positions  5,  6,  7,  were  reserved  for 
temperature  channel  selection.   For  example,  outputting  a 
binary  "16"  selected  channel  1,  outputting  a  binary  "64" 
selected  channel  4,  and  so  on.   Since,  by  design,  a  binary 
"0"  was  meant  to  reset  any  previous  binary  output  commands 
in  the  data  aquisition  system,  and  since  the  AD2  03  6  DPM 
calls  its  first  channel:  channel  "0",  any  time  a  binary  "0" 
is  output,  the  temperature  DPM  will  select  channel  "0" 
whether  or  not  the  computer  is  intending  to  read  the  AD2  03  6 
DPM. 
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IV.   DATA  AOUISITION 

A.   AQUIRING  DIGITAL  PANEL  METER  INFORMATION 

Figure  11  shows  the  general  procedure  for  aquiring  data 
from  one  of  the  DPM's.   Figure  12  illustrates  the  data 
conversion  procedure  for  obtaining  engineering  units. 
Adapter  "0"  was  reserved  to  control  the  scanivalve, 
turntable  position,  and  temperature  channel.   Adapter  "1" 
was  reserved  to  control  multiplexer  channel  selection. 
Control  sequence  for  reading  a  DPM  goes  as  follows:  (1) 
select  adapter  "1";  (2)  through  the  digital  output  device, 
output  a  binary  number  corresponding  to  the  multiplexer 
channel  handling  the  DPM  of  choice;  (3)  select  adapter  "0"; 
(4)  through  the  digital  input  device  (information  through 
the  multiplexer) ,  test  for  a  "data  ready"  signal  from  the 
DPM;  (5)  when  the  "data  ready"  signal  is  true,  read  the  DPM. 
The  DPM  will  have  ready  at  its  parallel  BCD  outputs,  the 
information  that  is  currently  displayed  on  the  face  of  the 
meter. 

Not  all  16  bits  of  the  digital  input  are  assigned  to 
numerical  information.   For  example,  the  highest  angle  of 
incidence  possible  would  be  360.0  degrees.   In  BCD  code, 
that  would  require  three  full  BCD  digits  and  one  "1/2"  BCD 
digit.   That  is,  the  most  significant  BCD  digit  would  only 
require  two  of  its  four  binary  bits  since  the  highest  number 
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required  of  that  BCD  digit  would  be  a  "3".   This  frees  two 
bits  of  information  that  could  be  used  for  other  purposes. 
Free  bits  are  devoted  to  such  uses  as  conveying  the  "data 
ready"  signal,  and  the  polarity  signal.   For  each  DPM,  both 
the  control  information  and  the  numerical  data  information 
inbound  to  the  computer  can  be  handled  through  one 
multiplexer  channel.   To  read  the  temperature  DPM,  add  step 
(3a) :  through  the  digital  outputs,  send  the  required  binary 
number  corresponding  to  the  temperature  channel  desired. 

All  of  the  DPM  information  already  exists  in  engineering 
units.   When  the  data  arrives  at  the  computer  through  the 
binary  input  port,  the  computer  does  not  "know"  that  it  is 
in  BCD.   After  the  data  is  read  from  the  binary  input  port 
into  a  variable  in  the  computer  program,  it  must  be 
converted  into  decimal  notation.   A  simple  algorithm  was 
written  to  accomplish  this. 

B.   AQUIRING  SCANIVALVE  INFORMATION 

Scanivalve  pressure  information  passes  through  the 
signal  conditioner/amplifier  and  then  directly  into  one  of 
the  analog  inputs  of  adapter  "0".   The  onboard  A/D  converter 
produces  an  integer  value  in  the  range  0-4  09  6.   This  number 
must  be  converted  into  engineering  units  in  software.   This 
is  accomplished  during  the  scanivalve  calibration  portion  of 
the  computer  program. 
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1.  Scanivalve  Calibration 

The  scanivalve  is  calibrated  by  applying  a  known 
pressure  source  to  the  scanivalve  transducer.   The  known 
pressure  source  was  provided  by  the  calibration  manometer. 
The  known  value  of  the  calibration  pressure  source  is  read 
from  the  calibration  pressure  DPM.   Both  transducers  are 
referenced  to  ambient  pressure.   When  the  computer  has 
obtained  the  above  information,  the  scanivalve  transducer 
can  be  calibrated  and  a  calibration  constant  calculated. 

The  scanivalve  is  commanded  to  two  ports,  one  for 
the  calibration  pressure,  and  the  second  for  reference,  or 
ambient  pressure.   The  information  returned  is  in  counts 
(calcount,  ref count,  from  the  A/D  converter) .   The  following 
equation  is  used  to  calculate  the  calibration  constant: 

calconstant  =  (calcount  -  refcount)  /  cal  pressure (PSI) . 

Calibration  is  complete  when  the  calibration  constant  is 
calculated.   The  units  of  the  calibration  constant  are 
counts  per  PSI. 

2 .  Scanivalve  Pressure  Information 

The  pressure  at  any  port  may  now  be  obtained  by 
commanding  the  scanivalve  to  that  port  and  obtaining  the  A/D 
count.   The  following  equation  is  then  used  to  convert 
pressure (count)  into  pressure (PSI) : 
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pressure (PSI)  =  (pressure (count)  -  refcount)  /  calconst 

3 .   Coefficient  of  Pressure  Calculation 

The  coefficient  of  pressure  is  obtained  by  first 
calculating  the  tunnel  dynamic  pressure:  Q.   Q  is  obtained 
by  commanding  the  scanivalve  to  the  ports  containing  the 
tunnel  settling  chamber  static  pressure  (P2)  ^rid  the  free 
stream  static  pressure  (P^) .   Q  is  calculated  by  the 
following  [Ref.  l:pp.  Al-4]: 


Q  =  ((P2(PSI)  -  P3^(PSI))  /  (1  -  (A^/A^)^))     *  1.075 


where  A2  =  area  of  the  test  section,  A-j_  =  area  of  the 
settling  chamber,  and  1.075  =  the  tunnel  calibration  factor, 
The  tunnel  calibration  factor  was  obtained  by  referencing 
^ideal'  i^i^^sured  through  the  two  static  pressure  manifolds, 
to  Q-true'  it^^^sured  through  a  special  calibrated  pitot  tube 
[Ref.  l:pp.  Al-5] .   The  coefficient  of  pressure  (Cp) 
follows: 

Cp  =  (pressure (PSI)  -  pstat(PSI))  /  Q. 
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V.   DATA  ANALYSIS 

A.  COMPUTER  MATHEMATICAL  INTEGRATION 

The  data  reduction  portion  of  the  program  requires 
various  integrations  to  be  performed.   Since  the  pressure 
data  can  be  varying  widely,  the  common  trapezoidal  rule  was 
avoided.   The  following  method  was  developed  to  integrate 
the  pressure  coefficient:  (1)  A  clamped  cubic  spline 
algorithm  splines  the  data  [Ref.  2:pp.  122-123].   To  use  the 
clamped  cubic  spline,  derivatives  at  the  end  points  are 
calculated  first.   The  spline  algorithm  produces  a  third 
degree  polynomial  between  each  pair  of  data  points;  (2) 
Using  Simpson's  Rule  [Ref.  2:pp.  188],  each  polynomial  is 
integrated  over  the  interval  between  each  pair  of  data 
points;   (3)  The  separate  integral  subtotals  are  summed  to 
provide  a  final  value. 

B.  DATA  REDUCTION 

1.   Force  Coefficients 

Cp  is  calculated  for  each  port  on  the  wing  model. 
By  the  integration  procedure  outlined  in  the  previous 
section,  the  force  coefficient  in  the  Y-direction  (Cy) ,  and 
the  force  coefficient  in  the  X-direction  (C^^)  can  be 
calculated  [Ref.  l:pp.  B2-5] .   The  following  equations 
apply: 
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Cy  = 


(2  ^*^P1  "  ^Pu)  d(x/c) 


Cx  =   f  ^""^   (Cpf  -  Cp^)  d(y/c) 


where  Cp-j^  =  Cp  over  the  lower  surface;   Cp^  =  Cp  over  the 
upper  surface;   Cp^  =  Cp  over  the  front  surface;   Cp^^  =  Cp 
over  the  rear  surface;   (x/c)  =  x  coordinate  normalized  to 
the  chord  length;   (y/c)  =  y  coordinate  normalized  to  the 
wing  thickness. 

2 .   Moment  Coefficient 

In  a  similar  manner,  the  moment  coefficient  about 
the  leading  edge  (Cjyj-Lg)  can  be  calculated  using  the 
following  equation: 

y 

/  2/c  .  1 

^Mle  =  1    (Cpf-Cpj.)*(y/c)  d(y/c)  -  (   (Cp^-Cp^)  (x/c)  d(x/c) 


3 .   Coefficient  of  Lift  and  Drag 

The  angle  of  incidence  (alpha)  is  obtained  from  the 
angle  of  incidence  DPM.   The  coefficient  of  lift  (C-^)  and 
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coefficient  of  drag  (C^)  can  now  be  calculated  using  the 
following  equations  [Ref.  l:pp.  B2-7]; 

C-jL  =  Cy  *  cos  (alpha)  ~  ^x  *   sin  (alpha) 
C^  =  Cj^  *  cos  (alpha)  -  Cy  *  sin  (alpha) 
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VI.   RESULTS 

A.  COEFFICIENT  OF  PRESSURE 

Figure  13  shows  a  sample  output  of  a  data  run  taken  at 
10  degrees  angle  of  incidence  and  Reynolds  number  (Re)  equal 
to  7.5  X  10^.   The  data  reveals  a  pressure  leak  in  the  model 
connections  at  about  (x/c)  =  .72.   Otherwise,  the  pressure 
distribution  conforms  to  expectations. 

B.  LIFT  COEFFICIENT 

Figure  14  shows  the  results  of  a  series  of  data  runs 
performed  at  Re  =  7.5  X  10^,  from   -18°  to  +20°  angle  of 
incidence.   The  graph  compares  data  from  a  previous  lab 
setup,  but  under  approximately  similar  conditions.   Figure 
(14)  also  shows  the  more  established  NACA  data.   As  shown  in 
the  graph,  the  experimentally  obtained  C-j_  is  consistently 
below  that  of  the  NACA  results.   This  may  be  accounted  for 
by  several  reasons.   First,  the  Reynolds  number  of  the  NACA 
data  (at  2.9  X  10^)  is  much  higher  than  of  the  experimental 
data.   Lower  Reynolds  numbers  tend  to  lower  the  C-|_ .   Second, 
there  may  be  an  error  in  the  angle  of  incidence  calibration 
which  could  shift  the  experimental  results  either  left,  or 
right.   Third,  there  is  a  likely  possibility  that  the  wing 
models  are  not  exactly  similar.   Fourth,  the  leak  in  one  of 
the  pressure  lines  has  caused  an  erroneous  data  point  that 
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is  forcing  the  integral  of  the  Cp  over  the  upper  surface  to 
be  higher  than  anticipated.   The  reason  why  the  experimental 
data  is  further  from  the  NACA  data  than  the  previous  results 
is  more  puzzling.   Perhaps  the  differences  lie  in  the 
different  techniques  of  integration.   The  previous 
calculations  used  the  trapezoidal  rule.   Also,  the  previous 
data  was  measured  at  a  slightly  higher  Reynolds  number  (8.9  - 
X  10  ) .   This  could  account  for  the  upward  shift. 

C.   DRAG  COEFFICIENT 

Figure  15  shows  more  results  calculated  from  the  same 
data  runs  mentioned  in  the  previous  paragraph.   The 
experimental  results  appear  to  generally  agree  with  the  NACA 
data  but  some  discrepancies  are  noted.   The  following  are 
postulated  as  possible  reasons.   The  existence  of  negative 
drag,  located  about  C-^   =    .48,  suggests  the  overall 
difficulty  in  working  with  the  small  numbers  associated  with 
the  C^^  calculations.   Noise  present  in  the  system  destroys 
the  resolution  of  the  measurements  and  becomes  magnified  in 
the  integrations  used  in  calculating  C-^   and  Cy.   The 
breakdown  in  the  overall  shape  of  the  drag  bucket  about  C-j^  = 
.48  indicates  that  the  data  is  close  to  the  experimental 
error.   Lastly,  the  existence  of  the  leak  in  one  of  the 
pressure  lines  may  have  influenced  the  drag  toward  negative 
numbers  due  to  the  erroneous  higher  pressure  on  the  aft 
portion  of  the  upper  surface  of  the  wing  model.   The 
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experimental  data  appears  closer  to  the  NACA  data  than  the 
previous  results.   The  general  shape  has  improved  but  data 
in  the  C-i=  0.0  -  1.0  range  prohibits  a  more  positive 
endorsement.   One  final  comment:   the  NACA  curves  probably 
resulted  from  the  smoothing  over  of  many  data  runs.   These 
experimental  results  resulted  from  a  single  run  taken  to 
test  the  data  aquisition  system..  A  systematic  approach  to 
this  experiment,  coupled  with  experience  in  working  with  the 
various  experimental  variables  should  improve  the  overall 
quality  of  the  results. 
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VII.   CONCLUDING  REMARKS 

A  data  aquisition  system  has  been  designed  that 
incorporates  all  of  the  necessary  elements  to  complete  the 
system  from  transducer  to  final  output.   All  of  the  hardware 
requirements  including  transducers,  input  circuits,  signal 
conditioning,  transmission,  processing,  storage,  and  output 
have  been  assembled  and  integrated  to  ensure  faithful 
representation  of  the  original  measurable  quantity.   A 
computer  has  been  integrated  into  the  design,  including  the 
computer  interfaces  such  as  a  digital  multiplexer,  and 
analog  to  digital  converters.   A  computer  program  has  been 
written  to  bring  the  system  to  life,  and  to  perform  the 
later  stages  of  signal  processing  necessary  to  complete  the 
data  analysis  portion  of  the  system.   Hardware  and  software 
requirements  of  the  final  output  have  been  assembled  and 
include  online  display  and  hardcopy  output  of  tabular  data 
and  graphical  data.   The  lessons  learned  from  this  thesis 
project  will  be  of  value  to  forthcoming  improvements  to  the 
Low  Speed  Wind  Tunnel  Facility  at  the  Naval  Postgraduate 
School. 
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VIII.   RECOMMENDATIONS 

In  working  with  the  tunnel  apparatus,  one  deficiency 
became  clear:  the  turntable  positioning  system  is  inadequate 
for  precise  measurements.   First,  there  is  too  much  play  in 
its  chain  drive  assembly.   Angle  of  incidence  adjustments 
that  require  direction  reversals  become  blurred  in  the  chain 
slack.   The  relay  /  DC  motor  system  is  imprecise  and  must 
rely  on  software  timing  to  correctly  position  the  turntable. 
Also,  the  relays  were  not  designed  for  rapid  on-off 
operation,  allowing  small  increments  of  turntable 
adjustment.   A  precision  stepper  motor  and  new  drive  system 
are  recommended  changes. 

Efficiency  could  be  improved  greatly  with  the  addition 
of  a  computer  controlled  tunnel  speed  adjustment  device.   At 
present,  the  program  user  must  manually  adjust  the  tunnel  Q. 
Also,  the  flexibility  of  the  experiment  would  be  improved  by 
adding  options  to  vary  the  tunnel  Q  as  well  as  angle  of 
incidence  in  the  same  data  recording  session. 

The  Scanivalve  transducer  used  to  obtain  the 
experimental  results  has  had  known  problems  in  past 
experiments.   Another,  more  accurate  transducer  has  been 
ordered  that  should  improve  the  accuracy  of  the  data.   It  is 
recommended  that  when  the  newer  transducer  arrives,  the 
experiments  be  repeated. 
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The  data  aquisition  and  control  adapter  comes  equipped 
with  an  internal  12-bit  analog  to  digital  converter.   The 
manufacturer  of  the  data  aquisition  and  control  adapter 
(Cyborg  Corp.)  has  available  higher  resolution  A/D 
converters  that  can  easily  be  incorporated  into  the  present 
system.   Using  higher  resolution  A/D  converters  with  the  new 
scanivalve  transducer  may  improve  the  accuracy  of  the  drag 
measurements. 

The  data  aquisition  system  built  for  this  thesis  project 
began  with  a  given  assortment  of  equipment  that  reduced  the 
overall  flexibility  and  usefulness  of  the  system  from  the 
start.   The  author  has  since  learned  that  there  are  systems 
available  from  reputable  manufacturers  (such  as  Hewlett 
Packard,  Omega  Engineering,  Fluke)  that  incorporate  complete 
integrated  designs  from  transducer  to  final  output.   These 
systems  could  provide  the  basis  of  an  extremely  flexible  and 
capable  data  aquisition  system  with  room  to  spare  for  growth 
and  enhancement.   They  are  recommended. 
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APPENDIX  A 

'MAIN  PROGRAM 

'SCANPRO  IS  THE  PROGRAM  FOR  THE  SCANIVALVE  EXPERIMENT. 

KEY  OFF:  CLS 
GOSUB  INITIALIZE 

S$  =  "1":  SUBROUTINE  =  1:  ENDPRO$  =  "1":  ERRORR$  =  "NO" 
WHILE  ENDPRO$  <>  "0" 
GOSUB  GREETINGS 
IF  ENDPRO$  =  "0"  THEN  GOTO  1 
S$  =  "1" 
WHILE  S$  <>  "0" 

ON  SUBROUTINE  GOSUB  PROMPTAOA,_ 

CALIBRATE, _ 
SPEED, _ 
DISPLAYVAR,_ 
RECORD, _ 
DISPLAYDATA,_ 
FILEDATA,_ 
REDUCEPL0T1,_ 
REDUCEPL0T2,_ 
PLOTCP 
SUBROUTINE  =  SUBROUTINE  +  1 
IF  SUBROUTINE  =  11  THEN  S$  =  "0" 
WEND 

SUBROUTINE  =  1 
1   WEND 

COLOR  7,0:  CLS:  END 

INITIALIZE: 

'ROUTINE  INITIALIZE  INITIALIZES  CERTAIN  DATA  AQUISITION 

AND  CONTROL 
'ADAPTER  VARIABLES,  AND  OTHER  DATA  STORAGE  VARIABLES 

NEEDED. 

' INPUTS :  NONE 
'OUTPUTS:  NONE 

CTRL%  =  0:  STOR%  =  0:  STAT%  =  0:  MODE%  =  0:  HNDSHK%  =  0 
RATE  =  500 
XORMSK%  =  &H0 

DIM  DATAVAR%(50) ,  COUNTS (50),  PRESSURE (50) , 
WANTAOA (80),  WINDATA (3900) 
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DIMX(20),  Y(20),  L(20) ,  U(20),  Z(20),  B(20),  C(20), 

D(20),  G(3,4) 
DIM  XWINGU(20) ,  XWINGL(20),  YWINGU(20),  YWINGL(20), 

CPUPPER(20),  CMLE(80) 
DIM  CPLOWER(20) ,  CX(80),  CY(80),  CL1(80),  CD1(80), 

ALPHAl (80),  CMC4 1(80) 
DIMCMAC(80),  DATA3$(11),  DATA2(10),  DATAl(lO),  CP(38), 

XCOORD(38) 
DIM  PLOTDU1(20) ,  PLOTDL1(20),  PLOTDU2(20),  PLOTDL2(2  0) 
DIM  PLOTDU3(20) ,  PLOTDL3(20),  PL0TA0A(3),  PLOTCOUNT$ ( 3 ) 
'NEXT  LINE  CONTAINS  DATA  FOR  XWINGU ( ) 
2     DATA  0.0,  .0033,  .0083,  .0125,  .0225,  .0363,  .0535,  .1, 

.16,  .203,  .307,  .404,  -.501,  .550,  .6,  .7,  .8,  .9,  - 
1.0 
'NEXT  LINE  CONTAINS  DATA  FOR  XWINGL() 
DATA  0.0,  .0033,  .0083,  .01,  .0218,  .0317,  .0483,  .098, 

.153,  .2,  .304,  .402,  .5,  .55,  .60,  .7,  .802,  .902, 
1.0 
'NEXT  LINE  CONTAINS  DATA  FOR  YWINGU ( ) 
DATA  0.0,  .012,  .0168,  .0197,  .0262,  .0327,  .0392, 

.0538, 

.0678,  .0755,  .0880,  .0938,  .0944,  .0924,  .0882, 

.0717,  .0464,  .0196,  0.0 
'NEXT  LINE  CONTAINS  DATA  FOR  YWINGL() 
DATA  0.0,  .009,  .0128,  .0143,  .0193,  .0232,  .0277, 

.0386, 

.0468,  .0527,  .0610,  .0648,  .065,  .0632,  .0604, 

.0484,  .0304,  .0120,  0.0 
FOR  I  =  1  TO  19:  READ  XWINGU (I):  NEXT  I 
FOR  I  =  1  TO  19:  READ  XWINGL(I):  NEXT  I 
FOR  I  =  1  TO  19:  READ  YWINGU (I):  NEXT  I 
FOR  I  =  1  TO  19:  READ  YWINGL(I):  NEXT  I 
RESTORE  2 

XCOORD(l)  =  0:  READ  DUMMY 

FOR  I  =  2  TO  3  6  STEP  2:  READ  XCOORD(I):  NEXT  I 
READ  DUMMY 
FOR  I  =  3  TO  3  5  STEP  2:  READ  XCOORD(I):  NEXT  I 

'DEFINE  A  FUNCTION  THAT  WILL  BE  USED  FOR  CUBIC  SPLINE  / 

INTEGRATE 
DEF  FNSPLINE(I,JJ)=Y(I)+B(I)*(JJ-X(I) )+C(I)*(JJ 
-X(I) )^2+D(I)*(JJ-X(I) )^3 
RETURN 
*  ******************************************************************** 

ERRORR: 

'THIS  ERROR  ROUTINE  DETECTS  WHEN  THE  DATA  AQUISITION 
'AND  CONTROL  ADAPTER  RETURNS  AN  ERROR  CONDITION  THROUGH 
'THE  VARIABLE  STAT% .    IF  AN  ERROR  DOES  OCCUR, 
•THE  PROGRAM  IS  AUTOMATICALLY  ABORTED.  (NOT  USER 
'FRIENDLY!) 
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'INPUTS:  STAT% 
'OUTPUTS:  NONE 

IF  STAT%  <>  0  THEN 

SCREEN  0:  LOCATE  12,10 

PRINT  USING  "EXECUTION  ERROR  ###";  STAT% 
LOCATE  2  0,10 

PRINT  "CONSULT  THE  DATA  AQUISITION  AND  CONTROL 
PROGRAMMING  SUPPORT" 
LOCATE  21,10 

PRINT  "MANUAL  APPENDIX  FOR  INFORMATION  ON  THIS 
EXECUTION  ERROR." 

LOCATE  23,27:  PRINT  "ENTER  ANY  KEY  TO  CONTINUE:" 
ERRORR$  =  INPUT$(1):  S$  =  "0":  GOTO  1 
END  IF 
RETURN 

SVPORT : 

'S.V.  PORT  LOCATION  READS  THE  SCANIVALVE  BCD  PORT 
LOCATION 

'INPUTS:  NONE 
•  OUTPUTS :  PORT 

ADAPT%  =  1:  RAWVALUE%  =  1:  DEVICE%  =8 

GOSUB  BINOUT    'CALL  BOUS ,  ENABLE  MULTIPLEXER  CHANNEL 

# 1 :  S.V.  PORT 

ADAPT%  =  0:DEVICE%  =  8:  COUNT  =  1 

ANDMSK%  =  &H7F  'MASK  OUT  UPPER  9  BITS,  ALLOW  BITS  0-6 

GOSUB  BININ    'CALL  BINM 

BCD  =  DATAVAR%(1):  GOSUB  BCD    'CALL  BCD  TO  DECIMAL 

PORT  =  DECIMAL 
RETURN 

************************************************************ 
BCD: 

'BCD  TO  DECIMAL  CONVERTS  THE  BCD/BINARY  NUMBER  RECEIVED 

FROM  BINM  OR 
'THE  S.V.  PORT  DISPLAY  TO  DECIMAL  FOR  NORMAL  USE. 

'INPUTS:  BCD 

' OUTPUTS :  DECIMAL 

DECIMAL  =  0 

FOR  EXPONENT  =  12  TO  0  STEP  -4 

TRUNC  =  FIX(BCD  /  (2  ^  EXPONENT)) 

DECIMAL  =  DECIMAL  +  (10  ^  (EXPONENT  /  4)  *  TRUNC) 
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BCD  =  BCD  -  (TRUNC  *  (2  ^  EXPONENT)) 
NEXT  EXPONENT 
RETURN 
************************************************************ 

RESETT: 

'RESET  SENDS  A  BINARY  "0"  OUT  TO  NULLIFY  ANY  "ONE'S" 

PREVIOUSLY  SENT. 
'IT  IS  USED  TO  CLEAR  RELAY  COMMANDS  SUCH  AS  STEP,  HOME, 

AOA+ ,  AOA- . 

'INPUTS:  NONE 
'OUTPUTS:  NONE 

ADAPT%  =  0:  RAWVALUE%  =  0:  DEVICE%  =  8 
GOSUB  BINOUT 
RETURN 

HOME : 

'HOME  "HOMES"  THE  SCANIVALVE  TO  PORT  NUMBER  ONE. 

•INPUTS:  NONE 
'OUTPUTS:  NONE 

ADAPT%  =  0:  RAWVALUE%  =  2:  DEVICE%  =  8:  COUNT  =  2000 
GOSUB  BINOUT    'CALL  BOUS ,  SEND  THE  SCANIVALVE  TO  PORT 

48 
GOSUB  RESETT 

GOSUB  DELAY    'WAIT  TWO  SECONDS,  LET  THE  S.V.  CATCH  UP 
GOSUB  STEPP    'STEP  THE  SCANIVALVE  TO  PORT  1 
RETURN 

STEPP: 

'STEPP  "STEPS"  THE  SCANIVALVE  TO  THE  NEXT  PORT 

LOCATION  AND  UPDATES 
'THE  VARIABLE  "PORT". 

'INPUTS:  NONE 
' OUTPUTS :  PORT 

ADAPT%  =  O:  RAWVALUE%  =  1:  DEVICE%  =  8:  COUNT  =  10  0 

GOSUB  BINOUT    'STEP  SCANIVALVE  TO  NEXT  PORT 

COUNT  =  1:  GOSUB  DELAY:  COUNT  =  100 

GOSUB  RESETT 

GOSUB  DELAY    'WAIT  FOR  1/4  0  SECOND  (PREVENT  S.V. 

BURNOUT) 
GOSUB  SVPORT    'VERIFY  LOCATION  OF  SCANIVALVE 

RETURN 
************************************************************ 
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ANALOGIN: 

•READS,  AVERAGES  50  SAMPLES  FROM  THE  ANALOG  CHANNEL 
SPECIFIED  BY  CHANLO% 

' INPUTS :  CHANLO% 
'OUTPUTS:  ANADATA 

ADAPT%  =  O:  DEVICE%  =  9:  COUNT  =  1000 

IF  CHANLO%  <>  0  THEN  GOTO  3    'DELAY  ONLY  IF  READING 

THE  SCANI VALVE 
GOSUB  DELAY    'DELAY  FOR  1  SECOND  TO  LET  PRESSURE 
STABILIZE 
3   COUNT  =  50 

GOSUB  ANAIN    'READ  50  VALUES 

SUM  =  0 

FOR  I  =  1  TO  COUNT    'AVERAGE  THE  50  VALUES 

SUM  =  SUM  +  DATAVAR%(I) 
NEXT  I 

ANADATA  =  SUM  /  COUNT 
RETURN 

READPORT : 

'READPORT  POSITIONS  THE  SCANIVALVE  TO  A  SINGLE  PORT 

SPECIFIED  BY 
'"WANTPORT",  READS  THE  DATA,  AND  RETURNS  THE  DATA  IN 

"ANADATA"  AND  COUNTS 

'INPUTS:  WANTPORT 
'OUTPUTS:  ANADATA,  COUNTS ( ) 

GOSUB  SVPORT    'FIND  S.V.  CURRENT  POSITION 

IF  WANTPORT  <  PORT  THEN  GOSUB  HOME 

WHILE  WANTPORT  >  PORT 
GOSUB  STEPP 

WEND 

CHANLO%  =  0 

GOSUB  ANALOGIN 

COUNTS (PORT)  =  ANADATA 
RETURN 

************************************************************ 
BINOUT: 

'THIS  ROUTINE  OUTPUTS  THE  DECIMAL  VARIABLE  "RAWVALUE%" 

IN  BINARY 
'VIA  TERMINALS  BOO  -  B015 

•INPUTS:  ADAPT%,  DEVICE%,  RAWVALUE%  (HNDSHK%  AND  STAT% 

ARE  INITIALIZED  AT  PROGRAM  START) 
•OUTPUTS:  NONE 
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CALL  BOUS  (ADAPT%,  DEVICE%,  HNDSHK% ,  RAWVALUE%,  STAT%) 

GOSUB  ERRORR 
RETURN 
************************************************************ 

BININ: 

'THIS  ROUTINE  INPUTS  BINARY  SIGNALS  FROM  THE  TERMINALS 

BIO  -  BI15,  AND 
•RETURNS  THE  VALUE  IN  THE  ARRAY  DATAVAR% . 

'INPUTS:  ADAPT%,  DEVICE%,  ANDMSK% ,  COUNT.   (HNDSHK%, 

MODE%,  STOR%, 
•OUTPUTS:  DATAVAR% (  )  (NOTE:  OFTEN  ONLY  THE  FIRST  ITEM 

IN  THE  ARRAY  DATAVAR% ( )  IS  REQUESTED) 

CALL  BINM  (ADAPT%,  DEVICE%,  HNDSHK% ,  MODE%,  STOR% , 
ANDMSK%,  XORMSK%,  COUNT,  RATE,  DATAVAR%(1), 
STAT%) 

GOSUB  ERRORR 
RETURN 

DELAY: 

•THIS  ROUTINE  USES  THE  DATA  AQUISITION  AND  CONTROL 
ADAPTER  TO  DELAY  •'COUNT"  MILLISECONDS. 

•INPUTS:  ADAPT%,  COUNT.   (STAT%  IS  INITIALIZED  AT 

PROGRAM  START) 
•OUTPUTS:  NONE 

CALL  DELAY  (ADAPT%,  COUNT,  STAT%) 
GOSUB  ERRORR 
RETURN 

ANAIN : 

'THIS  ROUTINE  INPUTS  AN  ANALOG  VOLTAGE  VIA  TERMINALS 

A/DO  -  A/D3,  AND 
'RETURNS  A  COUNT  IN  THE  RANGE  (O  -  4096) .   A  READING  IS 

TAKEN  "COUNT" 
'NUMBER  OF  TIMES.  EACH  READING  IS  RETURNED  IN  THE 

ELEMENTS  OF  DATAVAR% . 
'CHANLO%  SELECTS  THE  ANALOG  CHANNEL  TO  BE  READ  (0  -  3). 

'INPUTS:  ADAPT%,  DEVICE%,  CHANLO%,  COUNT,  RATE. 

(CTRL%,  MODE%,  STOR% ,  RATE  -  ARE  INITIALIZED  AT 

PROGRAM  START) 
' OUTPUTS :  DATAVAR% (  ) 

CALL  AINM  (ADAPT%,  DEVICE%,  CHANLO% ,  CTRL% ,  MODE%, 
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STOR%,  COUNT,  RATE,  DATAVAR%(1),  STAT%) 

GOSUB  ERRORR 

RETURN 
************************************************************ 

READAOA: 

'READ  AOA  READS  THE  CURRENT  TURNTABLE  POSITION. 

'INPUTS:  NONE 
'OUTPUTS:  AOA 

ADAPT%  =  1:  RAWVALUE%  =  2:  DEVICE%  =  8 
GOSUB  BINOUT    'ENABLE  MULTIPLEXER  CHANNEL  #2:  AOA 
ADAPT%  =  0:  COUNT  =  1:  ANDMSK%  =  &H8000 
GOSUB  BININ    'DATA  READY? 
'NOTE:  IF  "DATA  READY?"  IS  TRUE,  THE  VALUE  RETURNED 

WILL  BE  &H8000 
WHILE  DATAVAR%(1)  =  0 

GOSUB  BININ 
WEND 

ANDMSK%  =  &H4000 
GOSUB  BININ    'READ  POLARITY 
'NOTE  IF  "POLARITY"  IS  TRUE,  THE  VALUE  RETURNED  WILL  BE 

16384  (&H4000) 
POLARITY  =  DATAVAR%(1) 
ANDMSK%  =  &H3FFF 

GOSUB  BININ    'READ  AOA  (FINALLY!) 
BCD  =  DATAVAR%(1) 
GOSUB  BCD 

AOA  =  DECIMAL  /  10 

IF  POLARITY  =  0  THEN  AOA  =  AOA  *  -1 
RETURN 
************************************************************* 

AOACONTROL: 

•AOA  CONTROL  POSITIONS  THE  TURNTABLE  TO  THE  LOCATION 

SPECIFIED  BY 
' "WANTDEG" 

•INPUTS:  WANTDEG 
'OUTPUTS:  NONE 

4   GOSUB  READAOA 

IF  (AOA  >  22)  OR  (AOA  <  -2  0)  THEN  GOTO  4 

WHILE  (CINT(ABS (WANTDEG  -  AOA)  *  10)  /  10)  >=  .3 

ABSDIFF  =  ABS (WANTDEG  -  AOA):  IF  ABSDIFF  >  38  THEN 

ABSDIFF  =  38 
ADAPT%  =  0:  DEVICE%  =  8 

IF  WANTDEG  >  AOA  THEN  RAWVALUE%=4 :  GOSUB  BINOUT 
IF  WANTDEG  <  AOA  THEN  RAWVALUE%=8 :  GOSUB  BINOUT 
IF  (CINT(ABSDIFF*10)/10)  >= . 4  THEN  TIMING=100  ELSE 


60 


TIMING=1 
COUNT=(ABSDIFF*TIMING)+190:IF  COUNT>4000  THEN 

COUNT=4  000 
GOSUB  DELAY  : GOSUB  RESETT 
COUNT  =  500 

GOSUB  DELAY  •  'ALLOW  TURNTABLE  TO  STABILIZE 
5      GOSUB  READAOA 

IF  (AOA  >  23)  OR  (AOA  <  -2  0)  THEN  GOTO  5 
WEND 
RETURN 

READCAL: 

•THIS  ROUTINE  READS  THE  CALIBRATION  PRESSURE  AND 

RETURNS  THE  VALUE  IN 
'"CALPRESS",  IN  PSIG. 

•INPUTS:  NONE 
'OUTPUTS:  CALPRESS 

ADAPT%  =  1:  RAWVALUE%  =  4:  DEVICE%  =  8 

GOSUB  BINOUT    'ENABLE  MULTIPLEXER  CHANNEL  #3,  CAL. 

PRESSURE 
ADAPT%  =  0  :  COUNT  =  1:  ANDMSK%  =  &H8  000 
GOSUB  BININ    •DATA  READY? 
'NOTE:  IF  "DATA  READY?"  IS  TRUE,  THE  VALUE  RETURNED 

WILL  BE  &H8000 
'IF  FALSE,  THE  VALUE  WILL  BE  0 . 
WHILE  DATAVAR%(1)  =  0 

GOSUB  BININ 
WEND 
ANDMSK%  =  &H7FFF:  GOSUB  BININ    'READ  CALIBRATION 

PRESSURE 
BCD  =  DATAVAR%(1):  GOSUB  BCD    'CALL  BCD-TO-DECIMAL 
CALPRESS  =  (DECIMAL  /  100)  *  .014223    'CONVERT  CM 
-WATER  TO  PSIG 
RETURN 
************************************************************ 

READTEMP: 

'THIS  ROUTINE  READS  ONE  OF  SIX  CHANNELS  OF  THE 
TEMPERATURE  METER. 

•INPUTS:  WANTEMP  (INTEGER  0-5) 
•OUTPUTS:  TEMPER 

'NOTE:  THE  ANALOG  DEVICES  DIGITAL  THERMOMETER  DEFAULTS 
TO  CHANNEL  0  UNLESS  ANOTHER  CHANNEL  IS  REQUESTED 
VIA  THE  COMPUTER,  OR  MANUALLY 

'SELECTED  ON  THE  FRONT  PANEL. 

IF  WANTEMP  >0  THEN  TEMPER=0:  GOTO  ENDREADTEMP  'INVALID 
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RESPONSE 
'IF  NO  THERMO-COUPLE  IS  CONNECTED  TO  THAT  CHANNEL,  SKIP 

ROUTINE. 
ADAPT%  =  1:  RAWVALUE%  =  8:  DEVICE%  =  8 
GOSUB  BINOUT    'ENABLE  MULTIPLEXER  CHANNEL  #5, 

TEMPERATURES 
RAWVALUE%  =  WANTEMP  *  16 
ADAPT %  =  0 

GOSUB  BINOUT    'SELECT  A  TEMPERATURE  CHANNEL 
COUNT  =  999:  GOSUB  DELAY    'WAIT  1/100  SEC,  LET  TEMP 

STABILIZE 
COUNT  =  1:  ANDMSK%  =  &H2  000 
GOSUB  BININ    'DATA  READY?  ■ 
'NOTE:  THE  "DATA  READY?"  SIGNAL  ON  THIS  METER  GOES 

"HIGH"  FOR  5  MSEC  EVERY 
'CONVERSION  CYCLE  BUT  REMAINS  "LOW"  OTHER  TIMES. 
WHILE  DATAVAR%(1)  =  &H2  000 

GOSUB  BININ 
WEND 

ANDMSK%  =  &H1FFF:  GOSUB  BININ    'READ  TEMPERATURE 
BCD  =  DATAVAR%(1):  GOSUB  BCD    'CALL  BCD-TO-DECIMAL 
TEMPER  =  DECIMAL  /  10 
ENDREADTEMP: 
RETURN 

READATMOS : 

'THIS  ROUTINE  READS  ATMOSPHERIC  PRESSURE  AND  RETURNS 

THE  VALUE  IN 
'"ATMOS",  IN  PSIA. 

•INPUTS:  NONE 
'OUTPUTS:  ATMOS 

ADAPT%  =  1:  RAWVALUE%  =  16:  DEVICE%  =  8 

GOSUB  BINOUT    'ENABLE  MULTIPLEXER  CHANNEL  #5,  ATMOS. 

PRESS . 
ADAPT%  =  0:  COUNT  =  1:  ANDMSK%  =  &H2  000 
GOSUB  BININ    'DATA  READY? 
'NOTE:  IF  "DATA  READY?"  IS  TRUE,  THE  VALUE  RETURNED 

WILL  BE  &H2000 
'IF  FALSE,  THE  VALUE  WILL  BE  0. 
WHILE  DATAVAR%(1)  =  0 

GOSUB  BININ 
WEND 
ANDMSK%  =  &H1FFF:  GOSUB  BININ    'READ  ATMOSPHERIC 

PRESSURE 
BCD  =  DATAVAR%(1):  GOSUB  BCD 
ATMOS  =  DECIMAL  /  100 
RETURN 
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CALIBRATE: 

•THIS  ROUTINE  DETERMINES  THE  CALIBRATION  CONSTANT  VIA  A 

REFERENCE 
•PRESSURE  FROM  PORT  1  AND  A  CALIBRATION  PRESSURE  FROM 

PORT  48. 

•INPUTS:  NONE 

•  OUTPUTS :  CALCONST 

COLOR  14,1:  CLS 

LOCATE  5,29:  PRINT  "SCANIVALVE  CALIBRATION" 

COLOR  10:  LOCATE  10,13      • 

PRINT  ••USING  THE  HAND-CRANK  ON  THE  CALIBRATION 

MANOMETER  NEXT^^ 
LOCATE  12,13 
PRINT  ••TO  THE  COMPUTER,  APPLY  A  PRESSURE  BETWEEN  40  - 

50 
LOCATE  14,13 
PRINT  ••CENTIMETERS  OF  WATER.   THE  COLUMN  OF  WATER  ON 

THE   "  ; 
COLOR  12 
PRINT  ••LEFT  '• 
COLOR  10 
LOCATE  16,13 
PRINT  •'SHOULD  RISE  HIGHER  THAN  THE  COLUMN  ON  THE 

RIGHT .  ••  ; 
COLOR  7:  LOCATE  23,28:  PRINT  •'ENTER  ANY  KEY  WHEN 

READY :  •• 

SS$  =  INPUT$(1):  LOCATE  2  3,28:  PRINT  •• 
II 

COLOR  29:  LOCATE  18,23 

PRINT  •'CALIBRATING  SCANIVALVE  TRANSDUCER  .  .  .  •• 
GOSUB  READCAL    •OBTAIN  CALIBRATION  PRESS  IN  PSIG 
•NOTE:  HOME  THE  SCANIVALVE  TO  PORT  4  8  TO  SAVE  TIME: 
ADAPT%  =  0:  RAWVALUE%  =  2:  DEVICE%  =  8:  COUNT  =  2  000 
GOSUB  BINOUT:  GOSUB  RESETT:  GOSUB  DELAY 
WANTPORT  =48:  GOSUB  READPORT    •OBTAIN  CALIBRATION 

COUNTS 
CAL  =  ANADATA 

WANTPORT  =  1:  GOSUB  READPORT    'OBTAIN  REFERENCE  COUNTS 
REF  =  ANADATA 

CALCONST  =  (CAL  -  REF)  /  CALPRESS 
RETURN 

DISPLAYVAR: 

'THIS  ROUTINE  DISPLAYS  TUNNEL  TEMP,  AMBIENT  TEMP, 

TUNNEL  TOTAL  PRESSURE, 
•TUNNEL  FREESTREAM  STATIC  PRESSURE,   TUNNEL  Q,  TEST 

SECTION  FREESTREAM 
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'VELOCITY,  BAROMETRIC  PRESSURE. 

•INPUTS: TEMPER,  CALCONST,  ATMOS ,  PRESSURE,  TUNTEMP, 

VELOCITY 
'OUTPUTS:  AMBTEMP 

COLOR  1,1:  CLS:  LOCATE  12,38:  COLOR  31:  PRINT  "WAIT.  . 
II 

• 

GOSUB  READATMOS    'OBTAIN  BAROMETRIC  PRESSURE 
WANTEMP  =  0:  GOSUB  READTEMP    'OBTAIN  AMBIENT  TEMP 
AMBTEMP  =  TEMPER 
GOTO  7 
CALVE L: 

***********  MINI  ROUTINE:  CALCULATE  VELOCITY  ********** 
WANTEMP  =  1:  GOSUB  READTEMP    'OBTAIN  TUNNEL  TEMP 
TUNTEMP  =  TEMPER 

WANTPORT  =  2:  GOSUB  READPORT    'OBTAIN  TOTAL  PRESSURE 
GOSUB  PRESSURE    'CONVERT  TO  PSIA 
PTOTAL  =  PRESSURE 
WANTPORT  =  3:  GOSUB  READPORT    'OBTAIN  FREESTREAM 

STATIC  PRESS 
GOSUB  PRESSURE    'CONVERT  TO  PSIA 
PSTAT  =  PRESSURE 
Q  =  (PTOTAL  -  PSTAT)  *  1.075        '1.075  =  TUNNEL 

CALIBRATION 
IF  Q  <  .001  THEN  VELOCITY  =  0:  GOTO  6 
'NOTE:  FACTOR  .99  02  DUE  TO  NONZERO  VELOCITY  IN  THE 

SETTLING  CHAMBER 
VELOCITY  =  SQR((2  *  Q  *  1715  *  (TUNTEMP  +  459.67))  / 

(PSTAT  *  .9902) ) 

6  RETURN 

7  COLOR  14,1:  CLS 

LOCATE  5,25:  PRINT  "EXPERIMENTAL  CONDITIONS:" 
COLOR  12 

LOCATE  8,15:  PRINT  "BAROMETRIC  PRESSURE:   ";  ATMOS; 
"PSIA" 

II  r 


LOCATE  10,15:  PRINT 

"PSIA" 

LOCATE  12,15:  PRINT 

PSTAT;  "PSIA" 

LOCATE  14,15:  PRINT 
iipii 

LOCATE  16,15:  PRINT 
iipii 

LOCATE  18,15:  PRINT 

LOCATE  2  0,15:  PRINT 


TUNNEL  TOTAL  PRESSURE:   ";  PTOTAL; 
•TUNNEL  FREE  STREAM  PRESSURE:   "; 
'AMBIENT  TEMPERATURE:   ";  AMBTEMP; 
'TUNNEL  TEMPERATURE:   ";  TUNTEMP; 


'TUNNEL  Q:   ";  Q;  "PS I" 
'FREE  STREAM  VELOCITY:   "; 
VELOCITY;  "FT/SEC" 
LOCATE  22,15:  PRINT  "CALIBRATION  CONSTANT:   ";  CALCONST 
COLOR  10:  LOCATE  2  4,20:  PRINT  "ENTER  1  TO  CONTINUE,  0 

TO  ABORT:  " 
S$  =  INPUT$(1) 
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RETURN 
************************************************************ 

PRESSURE: 

'THIS  ROUTINE  CONVERTS  THE  RAW  DATA  IN  COUNTS  (ANADATA) 

TO  PRESSURE  IN 
'  PSIA. 

'INPUTS:  ANADATA,  REF,  CALCONST,  ATMOS 
•OUTPUTS:  PRESSURE 

PRESSURE  =  ( (ANADATA  -  REF)  /  CALCONST)  +  ATMOS 

RETURN 
************************************************************ 

SPEED: 

'THIS  ROUTINE  ALLOWS  THE  USER  TO  OBSERVE  TUNNEL 

VELOCITY  AS  HE  MANUALLY 
'ADJUSTS  THE  PITCH  OF  THE  PROP  SYSTEM. 

'INPUTS:  TEMPER,  CALCONST,  ATMOS,  PRESSURE 
'OUTPUTS:  VELOCITY,  KNOTS  TO  SCREEN 

COLOR  3 1 , 4 : CLS : LOCATE  5,27:  PRINT  "ADJUST  TUNNEL  SPEED 

NOW." 
COLOR  14:  LOCATE  8,30:  PRINT  "TUNNEL  VELOCITY:  " 
GOSUB  READATMOS    'OBTAIN  BAROMETRIC  PRESSURE 
WANTPORT  =  3:  GOSUB  READPORT:    'OBTAIN  FREESTREAM 
STATIC  PRESS 

GOSUB  PRESSURE    'CONVERT  TO  PSIA 
PSTAT  =  PRESSURE 

8  WANTPORT  =  2:  GOSUB  READPORT     'READ  TOTAL  PRESSURE 
GOSUB  PRESSURE    'CONVERT  TO  PSIA 

PTOTAL  =  PRESSURE 

WANTEMP  =  1:  GOSUB  READTEMP    'OBTAIN  TUNNEL 

TEMPERATURE 
TUNTEMP  =  TEMPER 

Q  =  (PTOTAL  -  PSTAT)  *  1.07  5        '1.07  5  =  TUNNEL 

CALIBRATION 
IF  Q  <  .001  THEN  VELOCITY  =  0:  GOTO  9 
'NOTE:  .9902  DUE  TO  NONZERO  VELOCITY  IN  THE  SETTLING 

CHAMBER 
VELOCITY  =  SQR((2  *  Q  *  1715  *  (TEMPER  +  459.67))  / 

(PSTAT  *  .9902) ) 

9  KNOTS  =  .5924  *  VELOCITY 

COLOR  10:  LOCATE  11, 30: PRINT  USING 

"####.##"; VELOCITY;: PRINT  "  FT/SEC" 

COLOR  13:  LOCATE  13,30:  PRINT  USING  "####.##"; 
KNOTS;: PRINT  "  KNOTS" 
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COLOR  11:  LOCATE  2  3,10 

PRINT  "ENTER  ANY  SINGLE  KEY  WHEN  FINISHED  ADJUSTING 

TUNNEL  SPEED:" 

S$  =  INKEY$:  IF  S$  =  ""  THEN  GOTO  8 

LOCATE  5,27:  PRINT  "  " 

LOCATE  2  3,10 

PRINT  " 
II 

LOCATE  23,25 

PRINT  "ENTER  1  TO  CONTINUE,  0  TO  ABORT:  " 
S$  =  INPUT$(1) 
RETURN 

GREETINGS : 

'THIS  ROUTINE  GREETS  THE  USER  AND  PROVIDES  INITIAL 

ORIENTATION  AND 
' INSTRUCTIONS . 

•INPUTS:  NONE 
'OUTPUTS:  NONE 

CLS:  COLOR  31:  LOCATE  12,38:  PRINT  "WAIT.  .  ." 

WANTDEG=0:  GOSUB  AOACONTROL    'RESET  AOA  TO  0. 

COLOR  0,7:  CLS 

LOCATE    5,20:    PRINT    "AEROLAB  WIND  TUNN 

EL" 
LOCATE  7,25:  PRINT  "WING     EXPERIMENT" 
LOCATE  10,10:  COLOR  4 
PRINT  "WELCOME  TO  THE  AEROLAB  WIND  TUNNEL  FOR  THE  WING 

EXPERIMENT  USING" 
PRINT  "     AN  IBM  PC-AT  TO  CONTROL  ALL  OF  THE  DATA 

AQUISITION  ACTIVITIES.   YOU" 
PRINT  "     MUST  MANUALLY  START  THE  WIND  TUNNEL  AND 

MANUALLY  ADJUST  THE  TUNNEL" 
PRINT  "     SPEED  TO  OBTAIN  THE  DESIRED  WIND  VELOCITY. 

BUT  FIRST,  THE  COMPUTER" 
PRINT  "     WILL  CALIBRATE  THE  SCANIVALVE  TRANSDUCER." 
COLOR  1:  LOCATE  2  3,23 

PRINT  "ENTER  1  TO  CONTINUE,  0  TO  ABORT" 
ENDPRO$  =  INPUT$(1) 
RETURN 
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************************************************************ 

PROMPTAOA : 

•THIS  ROUTINE  PROMPTS  THE  USER  TO  INPUT  THE  DESIRED 
ANGLE  OF  INCIDENCE. 

•INPUTS:  NONE 
•OUTPUTS:  WANTAOAO 

CLOSE 

GOSUB  10:  GOTO  12 

10  COLOR  14,3:  CLS 

11  LOCATE    5,20:     PRINT"A   NGLE  OF  INCIDENC 

E" 
LOCATE  7,30:  PRINT  "SELECTION":  RETURN 

12  OPEN  "I",  1,  "AOADATA.DAT" 
1  =  1 

WHILE  NOT  EOF(l) 

INPUT  #1,  WANTAOAl 
WANTAOA(I)  =  WANTAOAl 
1  =  1  +  1 

WEND 

CLOSE  1 

COLOR  4:  LOCATE  10,5 

PRINT  "DO  YOU  WISH  TO  USE  THE  PREVIOUS  SELECTIONS  OF 
ANGLE  OF  INCIDENCE?  " 

PRINT  "  (ENTER  1 

FOR  YES,  0  FOR  NO) " 

LOCATE  5,75:  S$  =  INPUT$(1) 

IF  S$  =  "0"  THEN  COLOR  14:  CLS:  GOSUB  11:  GOTO  13 

CLS: 1=1: LOCATE  5, 25: COLOR  5 : PRINT"PREVIOUS  AOA 
SELECTIONS: ":LOCATE  7,1 

COLOR  4 

WHILE  WANTAOA(I)  <>  -99 

PRINT  WANTAOA  ( I )  ;  •• ,   ••  ; 
I  =  I  +  1  ' 

WEND 

COLOR  0:  LOCATE  2  3,10 

PRINT  "ENTER  1  TO  USE  THESE  AOA  SELECTIONS,  0  TO  ENTER 
NEW  SELECTIONS:  " 

S$  =  INPUT$(1) 

IF  8$  =  "1"  THEN  GOTO  14 

GOSUB  10 

13  COLOR  4:  LOCATE  10,20:  PRINT  "ENTER  DESIRED  ANGLE  OF 

INCIDENCE:  " 
LOCATE  11,30:  PRINT  "(LIMITS  ARE  +20  TO  -18)" 
COLOR  5:  LOCATE  12,30:  PRINT  "(ENTER  -99  TO  STOP)" 
COLOR  4:  LOCATE  15,20:  PRINT  "TOTAL  NUMBER  OF 

SELECTIONS  ENTERED:  " 
LOCATE  16,25:  PRINT  "(MAX  NUMBER  OF  CHOICES  =  77)" 
SELECTNUM  =  0:  SELECT  =  0 


67 


COLOR  1 

WHILE  SELECT  <>  -99 

LOCATE  15,60:  PRINT  "     ":  LOCATE  15,60:  PRINT 
SELECTNUM 

LOCATE  10,58:  PRINT  "        ": LOCATE  10,58:  INPUT 
"  ",  SELECT 

WANTAOA( SELECTNUM  +1)  =  SELECT 

SELECTNUM  =  SELECTNUM  +  1 

IF  SELECTNUM  >  76  THEN  SELECT  =  -99: 
WANTA0A(SELECTNUM+l)=-99 
WEND 

LOCATE  2  0,6 
PRINT"ENTER  1  TO  STORE  THESE  FOR  LATER  USE,  0  NOT  TO 

STORE  THEM:" 
LOCATE  21,6 
PRINT"  (NOTE:  OLD  FILE  WILL  BE 

DESTROYED.) " 
LOCATE  20,70:  S$  =  INPUT$(1) 
IF  S$  =  "0"  THEN  GOTO  14 
OPEN  "0",  1,  "AOADATA.DAT" 
1  =  1 
WHILE  WANTAOA(I)  <>  -99 

PRINT  #1,  USING  "###  ";  WANTAOA(I) 

1  =  1  +  1 
WEND 

PRINT  #1,  USING  "###  ";  WANTAOA(I) 
CLOSE  1 

14  COLOR  0:  LOCATE  2  3,10 

PRINT  " 
II 

LOCATE  23,23:  PRINT  "ENTER  1  TO  CONTINUE,  0  TO  ABORT:  " 

8$  =  INPUT$(1) 

IF  S$  =  "0"  THEN  GOTO  15 

CLS:  LOCATE  12,38:  COLOR  31:  PRINT  "WAIT.  .  ." 

WANTDEG=WANTA0A(1) :  GOSUB  AOACONTROL    'POSITION  TO 

FIRST  AOA 
COLOR  3 0 , 6 : CLS : LOCATE  10 , 30 : PRINT"START  TUNNEL 

NOW": LOCATE  12, 23: COLOR  12 
PRINT"SET  APPROXIMATE  Q  IN  CM  OF  WATER" : COLOR  0 : LOCATE 

23,23 
PRINT" ENTER  1  TO  CONTINUE,  0  TO  ABORT:  ":  S$=INPUT$(1) 

15  RETURN 
************************************************************ 

RECORD: 

'THIS  ROUTINE  CONTROLS  THE  ACTUAL  RUNNING  OF  THE 
EXPERIMENT  AND  DATA  RECORDING. 

•INPUTS:  WANTAOAO,  COUNTS () ,  PRESSURE,  Q,  VELOCITY, 

TUNTEMP,  AMBTEMP  PSTAT,  PTOTAL,  ATMOS 
•OUTPUTS:  WINDATAO,  ENDATA 
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K  =  0:  M  =  1  :  RERUN$  =  "FALSE" 

VELOCITYl  =  VELOCITY 

WHILE  WANTAOA(M)  <>  -99  AND  S$  <>  "0" 

CLS:  LOCATE  12,35:  PRINT  "RESETTING  AOA.  .  ." 

WANTDEG  =  WANTAOA(M) :  GOSUB  AOACONTROL 

CLS:  LOCATE  12,31:  PRINT  "TUNNEL  STABILIZING.  .  .' 

COUNT=5000:  GOSUB  DELAY    'WAIT  FOR  5  SECONDS, 

STABILIZE  TUNNEL 

GOSUB  CALVEL     'CALL  CALCULATE  VELOCITY 

IF  ABS (VELOCITYl  -  VELOCITY)  >  2  THEN 

SCREEN  0 

GOSUB  SPEED  'CALL  SET  TUNNEL  VEL. 

END  IF 

IF  S$  =  "0"  THEN  GOTO  17 

CALL  SCRNPLTl 

LOCATE  22,5:  PRINT  "AOA  =" ;WANTAOA (M) 

FOR  J  =  4  TO  3  9 

WANTPORT  =  J:  GOSUB  READPORT 

GOSUB  PRESSURE 

Q  =  (PTOTAL  -  PSTAT)  *  1.075 

IF  Q  <  .001  THEN  CP  =  0 :  GOTO  16 

CP  =  (PRESSURE  -  PSTAT)  /  Q 

16  XPOS  =  INT(XC00RD(J-3)  *  520)  +  100 

YPOS  =  INT(CP  *  30  +  180)  +  100 

PSET  (XPOS, YPOS):  CIRCLE  (XPOS , YPOS) , 2 

L  =  (K  *  50)  +  J:  WINDATA(L)  =  CP 

Sl$  =  INKEY$ 

IF  Sl$  <>  ""  THEN 

LOCATE  2  3,31:  PRINT  " 
II 

COLOR  15:  LOCATE  2  3,13 

PRINT  "ENTER  1  TO  CONTINUE,  5  TO 

REPEAT  SAME  RUN,  0  TO  ABORT" 

S$  =  INPUT$(1) 

IF  S$  =  "5"  THEN 

J=4  0 

ELSEIF  S$  =  "0"  THEN 

SCREEN  0 

J  =  40 

END  IF 

LOCATE  2  3,13 

PRINT  " 

II 

LOCATE  2  3,31:  PRINT  "ENTER  ANY  SINGLE 

KEY  TO  PAUSE" 
COLOR  14 
END  IF 

IF  ERRORR$  <>  "NO"  THEN  ERRORR$="NO" :  J  = 
40:  RERUN $  =  "TRUE" 
NEXT  J 
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IF  (S$  =  "5")  OR  (S$  =  "0")  THEN  GOTO  17 

DATA1(1)=Q:  DATAl (2) =VELOCITY: 

DATAl ( 3 ) =WANTAOA (M) :  DATAl ( 4 ) =TUNTEMP 

DATAl ( 5 ) =AMBTEMP :  DATAl ( 6 ) =CALCONST : 

DATA1(7)=CALPRESS:  1=1 

FOR  J  =  40  TO  46 

L  =  (K  *  50)  +  J 

WINDATA(L)  =  DATAl (I) 

1  =  1  +  1 

NEXT  J 

DATA2 ( 1 ) =ATMOS :  DATA2 ( 2 ) =PTOTAL :  DATA2 ( 3 ) =PSTAT : 

1  =  1 

FOR  J  =  1  TO  3 

L  =  (K  *  50)  +  J 

WINDATA(L)    =    DATA2(I):    1=1+1 

NEXT   J 

IF  RERUN$  =  "TRUE"  THEN  RERUN$  =  "FALSE":  GOTO 

17 

M=M+l:   K=K+1 

17    WEND 

ENDATA  =  K  *  50 

SCREEN  0:  COLOR  31,5:  CLS :  LOCATE  12,38:  PRINT  "WAIT.  . 
II 

WANTDEG  =  0:  GOSUB  AOACONTROL    'RESET  AOA  TO  0. 

COLOR  1,3:  CLS:  LOCATE  12,15 

PRINT  "FINISHED  RECORDING  DATA.   "; 

COLOR  20:  PRINT  "SHUT  DOWN  TUNNEL  NOW." 

LOCATE  2  3,23:  COLOR  0  :  PRINT  "ENTER  1  TO  CONTINUE,  0 

TO  ABORT:" 
S$  =  INPUT$(1) 
RETURN 

FILE DATA: 

'THIS  ROUTINE  FILES  THE  DATA  LOCATED  IN  THE  ARRAY 

WINDATAO  AND  STORES 
'IT  ON  A  SEQUENTIAL  FILE  UNDER  THE  NAME   "WINDATA.TXT" 

' INPUTS :  WINDATA ( ) ,  ENDATA 
' OUTPUTS :  NONE 

COLOR  14,6:  CLS:  LOCATE  5,3  0 

PRINT    "FILE  DATA" 

COLOR  12:  LOCATE  8,15 

PRINT  "ENTER  1  TO  SAVE  THE  RECORDED  DATA,  0  NOT  TO  SAVE 

IT:" 
LOCATE  9,15 
PRINT  "  (NOTE:  OLD  DATA  WILL  BE 

DELETED.) " 
LOCATE  8,70:  SAVEDATA$  =  INPUT$(1) 
IF  SAVEDATA$  =  "0"  THEN  GOTO  18 
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OPEN  "O",  1,  "WINDATA.TXT" 
FOR  I  =  1  TO  ENDATA 

PRINT  #1,  USING  "####.####  ";  WINDATA(I) 
NEXT  I 
CLOSE  #1 
18  RETURN 

DISPLAYDATA: 

'THIS  ROUTINE  DISPLAYS  THE  DATA  AS  RECORDED  IN  THE 
EXPERIMENT. 

•  INPUTS :  WINDATA ( )  ,  ENDATA 
'OUTPUTS:  INFORMATION  TO  THE  SCREEN 

COLOR    14,2:    CLS :    LOCATE    1,25:    PRINT    "RECORDED 

DATA" 
COLOR  4:  LOCATE  4,2 
LOCATE  4,10:  PRINT  "WING  CP" 
FOR  I  =  1  TO  18 

IF  I  <  10  THEN  Jl$  =  RIGHT$(STR$(I) ,1)  ELSE  Jl$  = 

RIGHT$(STR$(I) ,2) 
J2$  =  RIGHT$(STR$(I  +  18), 2):  ROW  =1+5 
LOCATE  R0W,2:PRINT"P(";J1$;") :": LOCATE 

'ROW,20:PRINT"P(";J2$;")  :" 
NEXT  I 
COLOR  4 :DATA3$(1)=" BAROMETRIC  PRESSURE (PSIA) :  ": 

DATA3$(2)="TOTAL  PRESSURE (PSIA) :  " 
DATA3$(3)="FREESTREAM  STATIC  PRESS(PSIA):  ": 

DATA3$(4)="TUNNEL  Q(PSIA) :  " 
DATA3$ (5) ="TUNNEL  VELOCITY (FT/SEC) :  ":  DATA3$ ( 6) ="ANGLE 

OF  INCIDENCE (DEG) :  " 
DATA3$ (7) ="TUNNEL  TEMPERATURE (F) :  ":  DATA3$ (8) ="AMBIENT 

TEMPERATURE (F) :  " 
DATA3$ (9)=" CALIBRATION  CONSTANT:  ": 

DATA3$ (10) ="CALIBRATION  PRESSURE (PSIA) :  " 
FOR  I  =  1  TO  10 

R0W=I+3:  LOCATE  ROW, 40:  PRINT  DATA3$(I) 
NEXT  I 

COLOR  0:  LOCATE  2  0,50:  PRINT "ENTER  1  TO  CONTINUE" 
LOCATE  21,50  :  PRINT"       5  TO  STOP" 

LOCATE  22,50  :  PRINT"       0  TO  ABORT" 

K  =  ENDATA  /  50 
FOR  I  =  1  TO  K 

IF  I=K  THEN  COLOR  0 : LOCATE  22, 50: PRINT" 

": LOCATE  21,50:  PRINT"       0  TO 
ABORT": COLOR  1 
COLOR  11:  LOCATE  3,34:  PRINT"RUN: " ;  I:  COLOR  1 
FOR  J  =  4  TO  21 

L  =  (((I-l)  *  50)  +  J):  ROW  =  J  +  2 

LOCATE  ROW, 8: PRINT  WINDATA (L):  LOCATE  ROW, 26: 
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PRINT  WINDATA(L+18) 
NEXT  J 
FOR  J  =  1  TO  3 

L  =  ((I  -  1)  *  50)  +  J:  ROW  =  J  +  3 
LOCATE  ROW, 70:  PRINT  USING  "###.###"; 
WINDATA(L) 
NEXT  J 
FOR  J  =  40  TO  46 

L  =  ((I  -  1)  *  50)  +  J:  ROW  =  J  -  33 
LOCATE  ROW, 70:  PRINT  USING  "####.###"; 
WINDATA(L) 
NEXT  J 

LOCATE  22,70:  S$  =  INPUT$(1) 
IF  S$  <>  "1"  THEN  I  =  K 
NEXT  I 
RETURN 

STARTSPLINE: 

'TO  USE  A  CLAMPED  CUBIC  SPLINE,  DERIVATIVES  AT  THE  END 

POINTS  MUST 
'BE  EVALUATED: 
'ESTABLISH  A  MATRIX:  G,  THAT  WILL  BE  USED  TO  SOLVE  FOR 

THE 
'EQUATION  OF  A  QUADRATIC  AT  EACH  END  POINT. 
'CI  +  C2X  +  C3X^2  =  Y(X),  WHERE  X,  Y  ARE  KNOWN,  CI, 2, 3 

ARE  UNKNOWN. 

' INPUTS :  X ( ) ,  Y ( ) ,  N 
'OUTPUTS:  FPl,  FP2 


START  =  1:  FINISH  =  3:  GOSUB  19 

MATRIX 
GOSUB  2  0 

CONSTANTS 
FPl  =  2  *  C3  *  X(l)  +  C2 

DERIVATIVE 


'FORM  FIRST 
•SOLVE  FOR 
'FIRST 


START  =  N-2:  FINISH  =  N:  GOSUB  19 

MATRIX 
GOSUB  2  0 

CONSTANTS 
FP2  =  2  *  C3  *  X(N)  +  C2 

DERIVATIVE 
RETURN    'NOW  EXECUTE  SPLINE! 


'FORM  SECOND 
'SOLVE  FOR 
'SECOND 


19  'FORM  MATRIX  TO  SOLVE  FOR  COEFFICIENTS  OF  A  QUADRATIC: 

L  =  1 

FOR  I  =  START  TO  FINISH 

G(L,1)  =  1:  G(L,2)  =  X(I):  G(L,3)  =  X(I)^2:  G(L,4) 
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=  Y(I) 


L  =  L  +  1 
NEXT  I 
RETURN 


2  0  'CRAMERS  RULE  TO  SOLVE  FOR  THE  QUADRATIC  CONSTANTS: 


Dl  =  G(l,l) 
D2  =  G(l,2) 
D3  =  G(l,3) 

DET  =  Dl  + 


* 

* 

D2 
* 
* 
* 

D2 
* 
* 
* 

D2 


(G(2,2) 

(G(2,3) 

(G(2,l) 

+  D3 

(G(2,4) 

(G(2,3) 

(G(2,l) 

+  D3)  / 

(G(2,2) 

(G(2,4) 

(G(2,l) 

+  D3)  / 


G(3,3) 
G(3,l) 
G(3,2) 


G(2,3) 
G(2,l) 
G(2,2) 


G(3,2) 
G(3,3) 
G(3,l) 


Dl  =  G(l,l)  *  (G(2,4)  *  G(3,3) 

D2  =  G(l,4)  *  (G(2,3)  *  0(3,1) 

D3  =  G(l,3)  *  (G(2,l)  *  G(3,4) 

C2  =  (Dl  +  D2  +  D3)  /  DET 

Dl  =  G(l,l)  *  (G(2,2)  *  G(3,4) 

D2  =  G(l,2)  *  (G(2,4)  *  G(3,l) 

D3  =  G(l,4)  *  (G(2,l)  *  G(3,2) 

C3  =  (Dl  +  D2  +  D3)  /  DET 
RETURN 

EXECUTESPLINE: 


G(2,3) 
G(2,l) 
G(2,4) 

G(2,4) 
G(2,l) 
G(2,2) 


G(3,4) 
G(3,3) 
G(3,l) 

G(3,2) 
G(3,4) 
G(3,l) 


'THIS  ROUTINE  TAKES  THE  DERIVATIVES  FROM  "START  CUBIC 

SPLINE",  X() ,  AND 
•Y(),  AND  COMPUTES  A  SPLINE  FOR  THAT  INTERVAL 

'INPUTS:  FPl,  FP2,  X(),  Y() 

'OUTPUTS:  B(),  C(),  D()  (CONSTANTS  FOR  FNSPLINE()  ) 

ALPHAl  =  3  *  (Y(2)  -  Y(l))  /  (X(2)  -  X(l))  -  3  *  FPl 
ALPHAN  =  3  *  FP2  -  3  *  (Y(N)  -  Y(N-l))  /  (X(N)  -  X(N 

-1)) 

DEF  FNALPHA(I)=3*( ( Y (I+l) * (X (I) -X (I-l) ) ) - ( Y (I) * (X ( I+l) 
-X(I-1)))+(Y(I-1)*(X(I+1)-X(I))))/((X(I+1) 
-X(I))*(X(I)-X(I-1))) 
L(l)  =  2  *  (X(2)  -  X(l)) 
U(l)  =  .5 

Z(l)  =  ALPHAl  /  (2  *  (X(2)  -  X(l))) 
B(l)  =  FPl 
FOR  1=2  TO  N-1 

L(I)  =  2  *  (X(I+1)  -  X(I-l))  -  (X(I)  -  X(I-l))  * 

U(I-l) 
U(I)  =  (X(I+1)  -  X(I))  /  L(I) 

Z(I)  =  (FNALPHA(I)  -  (X(I)  -  X(I-l))  *  Z(I-l))  / 
L(I) 
NEXT  I 

L(N)  =  (X(N)  -  X(N-l))  *  (2  -  U(N-l)) 
Z(N)  =  (ALPHAN  -  (X(N)  -  X(N-l))  *  Z(N-l))  /  L(N) 
C(N)  =  Z(N) 
FOR  L=N-1  TO  1  STEP  -1 
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C(L)  =  Z(L)  -  U(L)  *  C(L+1) 
B(L)=((Y(L+1)-Y(L))/(X(L+1)-X(L)))-((X(L+1) 

-X(L) )*(C(L+1)+2*C(L) )/3) 
D(L)  =  (C(L+1)  -  C(L))  /  (3  *  (X(L+1)  -  X(L))) 
NEXT  L 
RETURN 

INTEGRATE : 

'THIS  ROUTINE  WILL  INTEGRATE  THE  DATA  GIVEN  TO  IT  IN 

THE  FORM  OF  X,Y 
•COORDINATES  BY  CREATING  A  SERIES  OF  THIRD  DEGREE 

POLYNOMIALS  (VIA  THE 
•SPLINE  ROUTINE)  AND  THEN  INTEGRATING  EACH  POLYNOMIAL 

OVER   EACH  INTERVAL 
'AND  LASTLY  SUMMING  EACH  SUBTOTAL  FOR  A  FINAL  ANSWER. 

'INPUTS:  X(),  Y(),  N;  USES  THE  FUNCTION:  FNSPLINE ( ) 
' OUTPUTS :  INTEGRAL 

GOSUB  STARTSPLINE 
GOSUB  EXECUTESPLINE 

'USE  SIMPSON'S  RULE  (EXACT  FOR  THIRD  DEGREE 
POLYNOMIALS) 

INTEGRAL  =  0:  SUM  =0 

FOR  I  =  1  TO  N-1 

JJ  =  X(I):  HH  =  (X(I  +  1)  -  X(I))  /  2 

SUM  =  FNSPLINE ( I, JJ) 

JJ  =  X(I)  +  HH:  SUM  =  SUM  +  FNSPLINE ( I, JJ)  *  4 

JJ  =  X(I  +  1):  SUM  =  SUM  +  FNSPLINE ( I, JJ) 

INTEGRAL  =  INTEGRAL  +  ( (HH  /  3)  *  SUM) 

NEXT  I 
RETURN 
************************************************************ 

CALCY : 

'THIS  ROUTINE  CALCULATES  THE  FORCE  COEFFICIENT   IN  THE 
"Y"  DIRECTION 

'INPUTS:  CPUPPERO,  CPLOWER ( ) ,  XWINGU(),  XWINGL() 
'OUTPUTS:  CY,  [  X(),  Y(),  N,  FOR  "INTEGRATE"  ] 

FOR  I  =  1  TO  19 

X(I)  =  XWINGU(I) 

Y(I)  =  CRUPPER (I) 
NEXT  I 
N  =  19 

GOSUB  INTEGRATE 
CY  =  INTEGRAL 
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FOR  I  =  1  TO  19 

X(I)  =  XWINGL(I) 
Y(I)  =  CPLOWER(I) 
NEXT  I 
N  =  19 

GOSUB  INTEGRATE 
CY  =  INTEGRAL  -  CY    ' CY  =  INTEGRAL (CPLOWER  -  CRUPPER) 

dX/C 
CY(J+1)  =  CY 
RETURN 
*  ******************************************************************** 

REDUCEPLOTl: 

•THIS  ROUTINE  CONTROLS  THE  DATA  REDUCTION  AND  PLOTTING 
DECISIONS. 

•INPUTS:  WINDATAO,  XWINGU ( ) ,  XWINGL(),  YWINGU ( ) , 

YWINGL ( ) ,  ENDATA 
•  OUTPUTS :  CRUPPER ( )  ,  CPLOWER ( ) 

•REDUCE  DATA  AND  PLOT  FOR  EACH  AOA  SEPARATELY: 

K  =  ENDATA  /  50 

COLOR  4,7:  CLS :  LOCATE  8,2  0 

PRINT  "ENTER  1  TO  PLOT  CP,  0  NOT  TO  PLOT  CP:  " 
PLOTFLAG$  =  INPUT$(1):  IF  PLOTFLAG$="0"  THEN  GOTO  21 
COLOR  1:  LOCATE  10,20:  LIMIT  =  K:  IF  LIMIT>3  THEN 

LIMIT=3 
PRINT  "YOU  MAY  PRINT  UP  TO";  LIMIT;  "PLOT(S)  OF  CP.": 

LOCATE  12,20 
COLOR  2:  PRINT  "ENTER  THE  NUMBER  OF  PLOTS  REQUESTED:" 
S$=INPUT$(1) :  PLOTNUM=  VAL(S$) : 
FOR  I  =  1  TO  PLOTNUM 

LOCATE  14,13:  COLOR  5 

PRINT  "ENTER  THE  AOA  CORRESPONDING  TO  THE  CP 

PLOT  REQUESTED:" 
LOCATE  14,70:  INPUT  "  ",  PLOTAOA(I) 
LOCATE  14,68:  PRINT  "        " 
NEXT  I 
PL0TC0UNT$(1)  =  "NO":  PLOTCOUNT$ (2 )  =  "NO": 

PL0TC0UNT$(3)  =  "NO" 
PLOTCOUNT  =  1 
COLOR  0:  LOCATE  23,23:  PRINT  "ENTER  1  TO  CONTINUE,  0 

TO  ABORT:" 
S$  =  INPUT$(1):  IF  S$  =  "0"  THEN  RETURN 

21      COLOR  14,1:  CLS : LOCATE  10,33:  PRINT  "REDUCING  DATA.  . 

II 

• 

FOR  J  =  0  TO  (K  -  1) 

•SEPARATE  OUT  UPPER  SURFACE  DATA  FROM  LOWER 
SURFACE  DATA: 
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CPUPPER(l)  =  WINDATA((J  *  50)  +  4) 

CPLOWER(l)  =  CPUPPER(l) 

CPUPPER(19)  =  WINDATA((J  *  50)  +  39) 

CPL0WER(19)  =  CPUPPER(19) 

FOR  I  =  2  TO  18 

CPUPPER(I)  =  WINDATA((J  *  50)  +  (2  *  I)  +  1) 
CPLOWER(I)  =  WINDATA((J  *  50)  +  (2  *  I)  +  2) 

NEXT  I 

'IF  PLOTTING  IS  DESIRED,  SET  UP  DATA  FILES 

IF  PLOTFLAG$  =  "0"  THEN  GOTO  2  4 
FOR  I  =  1  TO  PLOTNUM 

IF  WANTA0A(J+1)  =  PLOTAOA(I)  THEN 

PLOTCOUNT$(PLOTCOUNT)  =  "YES" 
PLOTCOUNT  =  PLOTCOUNT  +  1 

END  IF 
NEXT  I 

IF  PL0TC0UNT$(1)  <>  "YES"  THEN  GOTO  2  2 
FOR  I  =  1  TO  19 

PLOTDUl(I)  =  CRUPPER (I) 

PLOTDLl(I)  =  CPLOWER(I) 
NEXT  I 

22  IF  PL0TC0UNT$(2)  <>  "YES"  THEN  GOTO  2  3 

FOR  I  =  1  TO  19 

PL0TDU2(I)  =  CPUPPER(I) 

PL0TDL2(I)  =  CPLOWER(I) 
NEXT  I 

23  IF  PL0TC0UNT$(3)  <>  "YES"  THEN  GOTO  24 

FOR  I  =  1  TO  19 

PL0TDU3(I)  =  CPUPPER(I) 

PL0TDL3(I)  =  CPLOWER(I) 
NEXT  I 

24  COLOR  12,1:  LOCATE  12, 38: PRINT  "RUN:  ";  J+1 

GOSUB  CALCY 

GOSUB  CALCX 

ALPHA  =  (WINDATA( (J*50)  +42)  *  3.1415926#)  /  180 

CL  =  (CY  *  COS (ALPHA))  -  (CX  *  SIN (ALPHA)) 

CD  =  (CX  *  COS (ALPHA))  +  (CY  *  SIN (ALPHA)) 

GOSUB  CALCMLE 

GOSUB  CALBLOCK 
NEXT  J 
RETURN 
************************************************************ 

CALCX: 

'THIS  ROUTINE  CALCULATES  THE  FORCE  COEFFICIENT  IN  THE 
"X"  DIRECTION. 

'INPUTS:  CPUPPERO,  CPLOWER(),  YWINGU(),  YWINGL(), 
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-INTEGRAL 
•OUTPUTS:  CX,  [  X(),  Y(),  N,  FOR  "INTEGRATE"] 

FOR  I  =  1  TO  13  'UPPER  FRONT 

X(I)  =  YWINGU(I):  Y(I)  =  CPUPPER(I) 
NEXT  I 

N  =  13:  GOSUB  INTEGRATE 
CX  =  INTEGRAL 
FOR  I  =  1  TO  13  'LOWER  FRONT 

X(I)  =  YWINGL(I) :   Y(I)  =  CPLOWER(I) 
NEXT  I 

N  =  13:   GOSUB  INTEGRATE 
CX  =  CX  +  INTEGRAL 
FOR  I  =  1  TO  7  'UPPER  REAR 

X(I)  =  YWINGU(I  +  12):   Y(I)  =  CPUPPER(I  +  12) 
NEXT  I 

N  =  7:   GOSUB  INTEGRATE 
CX  =  CX  +  INTEGRAL 
FOR  I  =  1  TO  7  'LOWER  REAR 

X(I)  =  YWINGL(I  +12):   Y(I)  =  CPLOWER(I  +  12) 
NEXT  I 

N  =  7:  GOSUB  INTEGRATE 
CX  =  CX  +  INTEGRAL 
CX(J+1)  =  CX 
RETURN 

CALCMLE: 

'THIS  ROUTINE  CALCULATES  THE  COEFFICIENT  OF  MOMENT 
ABOUT  THE  LEADING  EDGE. 

'INPUTS:  CPUPPERO,  CPLOWER(),  XWINGU ( ) ,  XWINGL(), 

YWINGU(),  YWINGLO 
•OUTPUTS:  CMLE,  [  X(),  Y(),  N,  FOR  "INTEGRATE"  ] 

FOR  I  =  1  TO  19 

X(I)  =  XWINGU (I):   Y(I)  =  CPUPPER(I)  *  XWINGU (I) 
NEXT  I 

N  =  19:   GOSUB  INTEGRATE 
CMLE  =  INTEGRAL 
FOR  I  =  1  TO  19 

X(I)  =  XWINGL(I) :   Y(I)  =  CPLOWER(I)  *  XWINGL(I) 
NEXT  I 

N  =  19:   GOSUB  INTEGRATE 
CMLE  =  CMLE  -  INTEGRAL 
FOR  I  =  1  TO  13 

X(I)  =  YWINGU(I):   Y(I)  =  CRUPPER (I)  *  YWINGU(I) 
NEXT  I 

N  =  13:   GOSUB  INTEGRATE 
CMLE  =  CMLE  +  INTEGRAL 
FOR  I  =  1  TO  13 
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X(I)  =  YWINGL(I):   Y(I)  =  CPLOWER(I)  *  YWINGL(I) 
NEXT  I 

N  =  13:   GOSUB  INTEGRATE 
CMLE  =  CMLE  -  INTEGRAL 
FOR  I  =  1  TO  7 

X(I)  =  YWINGU(I+12) :   Y(I)  =  CPUPPER(I+12)  * 
YWINGU(I+12) 
NEXT  I 

N  =  7:   GOSUB  INTEGRATE 
CMLE  =  CMLE  +  INTEGRAL 
FOR  I  =  1  TO  7 

X(I)  =  YWINGL(I+12) :   Y(I)  =  CPL0WER(I+12)  * 
YWINGL(H-12.) 
NEXT  I 

N  =  7:   GOSUB  INTEGRATE 
CMLE  =  CMLE  -  INTEGRAL 
CMLE(J+1)  =  CMLE 
RETURN 

CALBLOCK: 

'THIS  ROUTINE  CALCULATES  SOLID  BLOCKING,  WAKE  BLOCKING, 

FOR  2-D  MODELS 
•AND  CORRECTS  THE  AERODYNAMIC  PARAMETERS:  CL,  CD,  Q, 

CMC/ 4. 

'INPUTS:  CL,  CD,  WINDATAO 

'OUTPUTS:  CLl,  CDl,  CLl () ,  CD1(),  ALPHA 1 ,  ALPHA1(), 
CMC41,  CMC41() 

EPSILON  =  .004981  +  (CD  /  (2  *  3.75)) 

'NOTE:  .004981  =  (Kl  *  MODEL  VOLUME)  /  (C^l.5)    [LAB 

MANUAL  PG  8-2] 
SIGMA  =  ((3.14159)^2  /  48)  *  (1  /  3.75)^2        '[LAB 

MANUAL  PG  8-4] 
CLl  =  CL  *  (1  -  SIGMA  -  (2  *  EPSILON)) 

' CL1=C0RRECTED  VALUE  OF  CL 
CLl (J  +  1)  =  CLl 
CDl  =  CD  *  ((1  -  (3  *  .004981)  -  (2  *  (CD  /  (2  * 

3.75))))) 
CDl (J  +  1)  =  CDl 
' CD1=C0RRECTED  VALUE  OF  CD 
Ql  =  WINDATA( (J+50)  +  40)  *  (1  +  (2  *  EPSILON))' 

Q1=C0RRECTED  VALUE  OF  Q 
CMC4  =  CMLE  +  (CL  /  4)       'CMC4=M0MENT 

COEFFICIENT  ABOUT  QUARTER  CHORD 
ALPHA  =  WINDATA( (J*50)  +  42) 
ALPHA 1  =  ALPHA  + 

( ( (57. 3*SIGMA)/ (2*3. 14159) ) * (CL+ (4*CMC4 ) ) ) 
ALPHA1(J+1)  =  ALPHAl 

'ALPHA1=C0RRECTED  VALUE  OF  AOA 
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CMC41  =  (CMC4  *  (1  -  (2  *  EPSILON) ) )  +  ((SIGMA  *  CLl)  / 

4) 
CMC41(J+1)  =  CMC41  'CMC41=C0RRECTED 

VALUE  OF  CMC4 
RETURN 

REDUCEPL0T2 : 

•THIS  ROUTINE  REDUCES  THE  THE  DATA  THAT  REQUIRES  INPUT 

FROM  EACH  SEPARATE 
'AOA  RUN;  HANDLES  PLOTTING  AND  HARDCOPY  DECISIONS  FOR 

THIS  DATA  ALSO. 

'INPUTS:  CMLEO,  CY  ( )  ,  K,  ALPHAl  ( )  ,  CDl  ( )  ,  CLl  ( )  , 

CMC41() 
'OUTPUTS:  XAC,  CMAC ( ) ;   DATA  TO  THE  PRINTER 

IF  K  =  1  THEN  GOTO  2  5    'NO  CALCULATIONS  IF  ONLY  ONE 

RUN 
'COMPUTE  THE  AERODYNAMIC  CENTER  (XAC)  BY  FINDING  THE 

SLOPE  OF 
'CMLE  VS  CY,  BY  USING  LEAST  SQUARES  FIT  [LAB  MANUAL  PG 

B2-10] : 

SUMl  =  0:  SUM2  =  0:  SUM3  =  0:  SUM4  =  0 
FOR  I  =  1  TO  K 

SUMl  =  SUMl  +  (CY(I)  *  CMLE (I)) 

SUM2  =  SUM2  +  CY(I) 

SUMS  =  SUM3  +  CMLE (I) 
SUM4  =  SUM4  +  (CY(I))^2 
NEXT  I 

DENOMINATOR  =  ( (K  *  SUM4 )  -  (SUM2^2)) 
IF  DENOMINATOR  =  0  THEN 

XAC  =  0 

ELSE  XAC  =  -(((K  *  SUMl)  -  (SUM2  *  SUM3))  / 
DENOMINATOR) 
END  IF 
FOR  I  =  1  TO  K 

CMAC (I)  =  CMLE (I)  +  (CLl (I)  *  XAC) 
NEXT  I 
25    COLOR  14,4:  CLS :  LOCATE  12,2 

PRINT  "ENTER  1  TO  OUTPUT  A  HARDCOPY  OF  CP,  CY ,  CX,  XAC, 

CMLE,  CL,  CD,  CMC/ 4,  CMAC," 
LOCATE  14,2 
PRINT  " 

ENTER  0  NOT  TO  OUTPUT : " 
HARDCOPY$  =  INPUT$(1) 

IF  HARDCOPY$  =  "0"  THEN  GOTO  2  6    'SKIP  HARDCOPY  OUTPUT 
FOR  I  =  1  TO  K 

LPRINT  "AOA:  "; ALPHAl (I) 

LPRINT  "PORT:",  "    CP:" 


79 


FOR  II  =  1  TO  3  6 

LPRINT  II,  WINDATA( (I-l)  *  50  +11  +  3) 
NEXT  II 

LPRINT:  LPRINT 
LPRINT  "CY=",  CY(I):  LPRINT  "CX=" ,  CX(I):  LPRINT 

"XAC=",  XAC 
LPRINT  "CMLE=",  CMLE(I):  LPRINT  "CL=",  CLl(I): 

LPRINT  "CD=",  GDI (I) 
LPRINT  "CMC/4=",  CMC41(I):  LPRINT  "CMAC=",  CMAC(I) 
LPRINT  CHR$(12) 
NEXT  I 
26    COLOR  7,0:  CLS 
RETURN 

PLOTCP: 

'THIS  ROUTINE  TAKES  THE  ARRAY  VARIABLES  SET  UP  WITH  THE 

X-Y  COORDINATES 
•AND  CALLS  THE  GRAPH  PROGRAM  VIA  THE  SHELL  STATEMENT 

•  INPUTS : 
' OUTPUTS : 

CLOSE 

IF  PLOTFLAG$  =  "0"  THEN  GOTO  27 
OPEN  "O",  #1,  "\GRAPH\PLOTDATU.DAT" 
FOR  I  =  1  TO  19 

PRINT  #1,  USING  "#.####^^^^    ";  XWINGU(I); 
PLOTDUl(I) 
NEXT  I 
CLOSE  #1 

OPEN  "0",  #1,  "\GRAPH\PLOTDATL.DAT" 
FOR  I  =  1  TO  19 

PRINT  #1,  USING  "#.####^^^^    "7  XWINGL(I); 
PLOTDLl(I) 
NEXT  I 
CLOSE  #1 
SHELL  "CPPLOT.BAT" 

IF  PLOTNUM  <  2  THEN  GOTO  27 

OPEN  "O",  #1,  "\GRAPH\PLOTDATU.DAT" 

FOR  I  =  1  TO  19 

PRINT  #1,  USING  "#.####----    ";  XWINGU(I); 
PL0TDU2(I) 
NEXT  I 
CLOSE  #1 

OPEN  "O",  #1,  "\GRAPH\PLOTDATL.DAT" 
FOR  I  =  1  TO  19 

PRINT  #1,  USING  "#.####^^^^    ";  XWINGL(I); 
PL0TDL2(I) 
NEXT  I 
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CLOSE  #1 

SHELL  "CPPLOT.BAT" 

IF  PLOTNUM  <  3  THEN  GOTO  27 

OPEN  "O",  #1,  "\GRAPH\PLOTDATU.DAT" 

FOR  I  =  1  TO  19 

PRINT  #1,  USING  "#.####^^^^    ";  XWINGU(I); 
PL0TDU3 (I) 
NEXT  I 
CLOSE  #1 

OPEN  "0",  #1,  "\GRAPH\PLOTDATL.DAT" 
FOR  I  =  1  TO  19 

PRINT  #1,  USING  "■#.####'-^^^    ";  XWINGL(I); 
PL0TDL3(I) 
NEXT  I 
CLOSE  #1 

SHELL  "CPPLOT.BAT" 
27    RETURN 

SUB  SCRNPLTl   STATIC 

GOSUB  DIMARR 
GOSUB  HIRES 
GOSUB  VERTICALGRID 
GOSUB  HORIZONTALGRID 
GOSUB  NUMBERY 
GOTO  ENDSCRNPLTl 
************************************************************ 

DIMARR: 

DIM  GRID(IO) 

DIM  XGRID(20) ,  YGRID(IO),  LCHAR$(65),  E(3) 
RETURN 

HIRES: 

SCREEN  9  :  COLOR  14,5:  CLS 

WINDOW  (0,0)  -  (639,349) 

LINE  (100,  100)  -  (620,  340),,  B 
RETURN 
************************************************************ 

VERTICALGRID: 

COLOR  10 
FOR  I  =  1  TO  7 

XGRID(I)  =  I  *  65  +  100 

LINE  (XGRID(I),  340)  -  (XGRID(I),  98),,,&HEEEE 
NEXT  I 

LOCATE  19,13:  PRINT  "0.0":  LOCATE  19,28:  PRINT 

"0.25" 
LOCATE  19,45:  PRINT  "0.5":  LOCATE  19,61:  PRINT 
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"0.75" 
LOCATE  19,77:  PRINT  "1.0" 

COLOR  12;  LOCATE  21,38:  PRINT  "WING  CHORD  (X/C) " 
COLOR  15:  LOCATE  2  3,31:  PRINT  "ENTER  ANY  SINGLE 

KEY  TO  PAUSE" 
RETURN 

HORIZONTALGRID: 


GRID (1) =-6: GRID (2) =-4: GRID (3)= 

-2 : GRID ( 4 ) =0 : GRID ( 5 ) =2 : GRIDNUM=5 
SCALEY  =  (GRID(GRIDNUM)  -  GRID(l))  /  240 
YCONSTANT  =  GRID(l)  /  SCALEY 
COLOR  10 
FOR  I  =  1  TO  GRIDNUM 

YGRID(I)  =  INT (GRID (I)  /  SCALEY  -  YCONSTANT)  +  100 
IF  GRID(I)  =  0  THEN  LINE  (98 , YGRID (I) ) 

-(620,YGRID(I) ) :  GOTO  28 
IF  (1=1)  OR  (I=GRIDNUM)  THEN  GOTO  2  8 
LINE  (98,  YGRID(I))  -  (620,  YGRID (I) ) , , , &HFOFO 
28    NEXT  I 

COLOR  12:  LOCATE  10,5:  PRINT  "CP" 
RETURN 

SIDESTRING: 

'ROUTINE  TO  PRINT  A  STRING  AT  A  SPECIFIC  SCREEN 
POSITION,  SIDEWAYS 


VARIABLES:  WORD$   - 
LONG 
C$ 
WORD$ 
ASCII   - 
X 
CHARACTER 

Y 
CHARACTER 
INPUTS:  WORD$,  X,  Y 
OUTPUTS:  NONE 


STRING  TO  BE  PRINTED 
LENGTH  OF  WORD$ 
INDIVIDUAL  CHARACTERS  OF 

ASCII  DECIMAL  VALUE  OF  C$ 
X  POSITION  OF  FIRST 

Y  POSITION  OF  FIRST 


'DEFINE  THE  SIDEWAYS  CHARACTER  SET: 
COLOR  14 

LCHAR$(1)  =  "BU7" 
'BLANK 
LCHAR$(9)  =  "BL2  L5  H2  BR9  Gl  BU3  BRl  BU2" 

'( 
LCHAR$(10)=  "BL9  BU2  E2  R5  F2  BU4" 

•) 

LCHAR$(14)=  "BUI  BL5  BUI  U3  BR5  BU2" 
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LCHAR$(15)=  "UO  BU2" 
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LCHAR$(17)= 

LCHAR$(18)= 

LCHAR$(19)= 

LCHAR$(20)= 

LCHAR$(21)= 

LCHAR$(22)= 

LCHAR$(23)= 

LCHAR$(24)= 

LCHAR$(25)= 

LCHAR$(2  6)= 

LCHAR$(34)= 

LCHAR$ (35)= 

LCHAR$ (36)= 

LCHAR$(37)= 

LCHAR$(38)= 

LCHAR$(39)= 

LCHAR$(40)= 

LCHAR$(41)  = 

LCHAR$ (42)= 

LCHAR$ (43): 

LCHAR$(44): 

LCHAR$(45) 

LCHAR$ (46) 

LCHAR$(47) 

LCHAR$(48) 


"BLl  L7  HI  U2  El  R7  Fl  ND2  BU3" 

0 

"BU2  NL9  BU2" 

1 

"BL8  LI  U4  R2  F4  R3  U4  BU2" 

2 

"BL8  HI  U2  El  R2  Fl  ND2  El  R3  Fl  D2  NGl 
BU5"    '3 

"BL4  L5  BU4  R5  ND3  R4  BU2" 

4 

"BL5  NU3  L4  U4  BR4  R5  Fl  D2  Gl  BU4  BRl  BU2" 

5 

"BLl  L7  HI  U2  BR4  ND3  El  R3  Fl  ND2  BU3" 

6 

"BL9  U4  R5  F4  BU6" 

7 

"BLl  L3  BLl  L3  HI  U2  El  R2  Fl  ND2  El  R3  Fl 
ND2  BU3"  '8 

"BL6  L2  HI  U2  El  R3  ND3  R5  BU2" 

9 

"L7  H2  E2  R2  ND3  R5  BU2" 

A 

"L9  U3  El  R3  ND3  R4  Fl  ND2  BU3" 

B 

"BLl  L7  HI  U2  El  BR6  Rl  Fl  ND2  BU3" 

C 

"L9  U3  El  R7  Fl  ND2  BU3" 

D 

"L9  NU4  BR4  NU3  BR5  U4  BU2" 

E 

"L9  NU4  BR4  NU3  BR5  BU6" 

F 

"BLl  L7  HI  U2  El  BR3  ND4  R4  Fl  ND2  BU3" 

G 

"L9  BU4  R4  ND3  R5  BU2" 

H 

"BL9  U2  NU2  R9  ND2  U2  BU2" 

I 

"BLl  L2  BU4  NL6  R2  Fl  ND2  BU3" 

J 

"L9  BU4  F4  E4  Rl  BU2" 

K 

"NL9  U4  BU2" 

L 

"L9  E2  NRl  H2  R9  BU2" 

M 

"L9  Ul  R3  Ul  R3  Ul  R3  Ul  NL9  BU2" 

N 

"BLl  L7  HI  U2  El  R7  Fl  ND2  BU3" 

0 
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LCHAR$(49)=  "L9  U3  El  R2  Fl  ND2  BUI  BR5  BU2" 

'P 
LCHAR$(50)=  "BLl  L7  HI  U2  El  R8  NG2  ND3  BU2" 

•Q 
LCHAR$(51)=  "L9  U3  El  R2  Fl  D2  E3  R2  BU2" 

'R 
LCHAR$(52)=  "BL8  HI  U2  El  BD4  R2  El  U2  El  R3  Fl  D2  Gl 

BU4  BRl  BU2"  'S 
LCHAR$(53)=  "BL9  U2  NR9  U2  BR9  BU2" 

•T 
LCHAR$(54)=  "BLl  NL8  BU4  NL8  BRl  BDl  ND2  BU3" 

'U 
LCHAR$(55)=  "BU3  Gl  LI  Gl  L6  BU4  R6  Fl  Rl  Fl  BU4" 

'V 
LCHAR$(56)=  "NL9  H2  NLl  E2  NL9  BU2" 

'W 
LCHAR$(57)=  "L3  H2  G2  L2  BU4  R2  Fl  BR2  El  R3  BU2" 

•X 
LCHAR$(58)=  "BU2  L5  G2  L2  BU4  R2  Fl  BR6  BU3" 

•Y 
LCHAR$(59)=  "BL9  LI  U4  R3  F4  R3  U4  NLl  BU2" 

•Z 

Y  =  -Y  +  349 

MOVE$  =  "BM"  +  STR$(X)  +  ","  +  STR$(Y) 
.  DRAW  "X"  +  VARPTR$(MOVE$) 
LONG  =  LEN  (WORD$) 

IF  NEGFLAG$  =  "Y"  THEN  START  =  1  ELSE  START  =  2 
FOR  J  =  START  TO  LONG 

C$  =  MID$(WORD$,  J,  1) 
ASCII  =  ASC(C$)  -  31 
DRAW  "X"  +  VARPTR$ (LCHAR$ (ASCII) ) 
NEXT  J 
RETURN 
•  ************************************************************* 

NUMBERY : 

•THIS  ROUTINE  WILL  TAKE  THE  VALUES  OF  THE  Y-AXIS 

GRIDLINES  AND 
'NUMBER  THE  Y-AXIS  WITH  SIDEWAYS  CHARACTERS. 

'INPUTS:  GRIDNUM,  YGRID() 
'OUTPUTS:  WORD$ 

=  96 
FOR  I  =  1  TO  GRIDNUM 
X  =  96 

Y  =  YGRID(I)  -  5:  IF  GRID(I)>=10  THEN  Y  =  Y  -  3 
WORD$  =  STR$(GRID(I) ) 

IF  GRID(I)<0  THEN  NEGFLAG$  =  "Y"  ELSE  NEGFLAG$  = 
"N" 
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IF  LEN(WORD$)  =  2  THEN  WORD$  =  STR$ (GRID (I) )  +  "." 

+  "0" 
GOSUB  SIDESTRING 
NEXT  I 
RETURN 
ENDSCRNPLTl: 
END  SUB 


85 


LIST  OF  REFERENCES 


1.  Laboratory  Manual  for  Low  Speed  Wind  Tunnel  Testing, 
Department  of  Aeronautics,  Naval  Postgraduate  School, 
October  1983. 

2.  Burden,  R.  L. ,  Faires,  J.  D. ,  and  Reynolds,  A.  C. , 
Numerical  Analysis.  Prindle,  Weber  &  Schmidt,  1978. 


86 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1.  Defense  Technical  Information  Center  2 
Cameron  Station 

Alexandria,  Virginia    22304-6145 

2.  Library,  Code  0142  2 
Naval  Postgraduate  School 

Monterey,  California    93943.-5002 

3.  Satyanarayana  Bodapati  10 
M.  S.  260-1 

NASA  Ames  Research  Center 
Moffet  Field,  CA    94035 

4.  M.  F.  Platzer,  Chairman  1 
Department  of  Aeronautics 

Code  67 

Naval  Postgraduate  School 

Monterey,  Ca    93943-5100 

5.  John  A.  DiMiceli  5 
12  65  Leahy  Road 

Monterey,  CA    93  94  0 

6.  Ted  Dunton  1 
Department  of  Aeronautics 

Code  67 

Naval  Postgraduate  School 

Monterey,  Ca    93943-5100 


87 


/v 


1 - 


WAv,:...i   '->«'••.  t^.i-....u^irg  SCHOOL 
MOIJTEBBY.  CALIFOEmA  93943  B002 


i 


2201G5 

^Thesis 
TiD57646     Dimiceli 
^cA  Computer-controlled 

data  acquisition  and 

analysis. 


8     FEB  90 


«r>T 


m 


H 


Thesis 
D57646 
c.l 


220185 


Dimiceli 

Computer-controlled 

data  acquisition  and 

analysis. 


